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Disclaimer 

Mention  of  trade  names  or  commercial  products  does  not  constitute  endorsement  or 
recommendation  for  use  by  the  government.  No  author  or  anyone  acting  on  behalf  of  the 
software  system  assumes  any  liabilif^  or  responsibility  for  damages  (including  consequen- 
tial) caused  by  reliance  on  the  materliis  presented,  or  from  the  use  of  any  information, 
apparatus,  or  methods  of  processing  disclosed  in  this  manual.  The  right  is  reserved  to 
make  changes  to  the  materials  contained  herein  w^ithout  notice. 
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Training  Guide  for 
Solid  IVIinerals 
Specialists 

In  response  to  a  growing  demand 
within  the  Bureau  of  Land  Management 
for  GIS  training  beyond  the  introductory  level, 
several  advanced  courses  have  been  developed. 
This  manual  is  to  be  used  in  training  BLM 
minerals  specialists  to  use  MOSS/MAPS  for 
sohd  minerals  apphcations.  While  this  text 
should  be  of  continued  value  to  GIS  users, 
it  is  not  meant  to  replace  more  complete 
documentation  and  user  manuals  already 
in  existence.  Rather,  this  is  a  training  manual 
and  a  quick  reference  document. 
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Chapter  Ones 

Systems  Operations 
A  Review 


In  this  chapter  you  will  review 

■  GIS  directory  structures 

■  MOSS-family  systems  interfaces 

■  PRIMOS  system  commands 

■  Disk  management 

■  Text  editor  (EMACS) 

■  Non-interactive  processing 

■  GIS  project  Management 

A  review  exercise  will  follow  this  chapter 
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Directory  Structure 


<MFD> 


Master  File  Directory  is  the  basic  unit  of  storage  on  the 
PRIME. 


<MFD>UFD>         User  File  Directory  is  a  sub-directory  of  the  Master  File 
Directory;  it  may  contain  many  files  and  other  sub- 
directories. 

Examples: 

1.  Pathname:  <MFD6>userl>MOSS 

MFD6 


User  1  User2,  User3 


MOSS  MOSS  MOSS 

2.  Pathname:  <2VIFD4>MOSSDATA>MINTRAIN 


MFD4 


MOSSDATA 


NEW.MEX    WOLSAG   MINTRAIN 


map  files       map  files       map  files 
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MOSS  Family  Systems  Interfaces 

Several  software  packages  (systems)  may  be  required  to  satisfy  BLM  GIS 
project  requirements.  For  example,  a  solid  mineral  application  of  a  GIS  may 
require  you  to  incorporate  the  Geological  Surveys  Digital  Line  Graph  (DLG) 
data  into  MOSS.  The  system  DLGS20  converts  DLG  standard  format  to 
optional  format,  and  then  the  system  DLG3  converts  the  DLG  optional 
format  to  a  data  file  that  can  be  imported  into  MOSS.  Therefore,  an  under- 
standing of  MOSS-family  system  interfaces  is  essential  to  GIS  project  plan- 
ning and  implementation. 

Figure  1.1  is  a  diagram  that  shows  the  various  systems  available  and  how 
they  interface  within  the  MOSS-family  software. 

Table  1.1  Hsts  the  system  names  and  their  functions. 
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Figure  1 .1   MOSS-Famiiy  Software  System  Interface 
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TABLE  1.1.  SYSTEM  INTERFACE  REFERENCE 


PltOGRAMNAME        FUNCTION 


LOCATION 


ADS.ADSTODLG 
ADS2MOSS 
ADS  GENERATE 

CNVADS 

CONTOUR 

CTOG 

DLG02ADS 

DLGS20 

DLG3 

EXPORT 
ADS.GGWL 
GRID 
IMPORT  (MAPS) 

IMPORT  (MOSS) 

MAPSTOIDIMS 

MOSS2ADS 

MOSS2LGS 

ADS.PARCEL 

PCCS2ADS 
PITOGGWL 

prroMoss 

POLYCELL 

RASTEREZE 

ADS.REV6T07 

SPSS 

XYSUBJECT 


Converts  ADS  Eles  to  DLG  optional  fonnat 

Converts  ADS  files  to  MOSS 

Produces  ADS  files  from  COGO  or  ASCH  files 
Also  allows  screen  digitizing. 

Converts  SAM  file  to  DAM  file 

Converts  cell  data  to  line  data 

Converts  MOSS  IMPORT/EXPORT  fonnat  to  DEM 

Converts  DLG  optional  format  to  ADS 

Converts  DLG  standard  format  to  optional 

Converts  DLG  optional  format  to 
MOSS  IMPORT/EXPORT  fonnat. 

Converts  MOSS  files  to  IMPORT/EXPORT  fonnat 

Converts  GGWL  files  to  ADS 

Converts  point  data  to  cell  data 

IDIMS  byte  or  word,  USGS  DEM, 
DMA  3  arc-second  DEM  to  MAPS. 

IMPORT/EXPORT  fonnat  to  MOSS 

MAPS  files  to  IDIMS  format 

MOSS  files  to  ADS 

MOSS  files  (type  1,2,3)  to  LGS  format 

ADS  landlines  and  ALMRS  status  or 
match/merge  files  to  ADS. 

Converts  PCCS  files  to  ADS 

EMACS-edited  Petroleum  Information  (PI) 
data  to  ADS.GGWL 

Petroleum  Information  data  to 
MOSS  IMPORT/EXPORT  fonnat. 

MOSS  (vector)  files  to  MAPS  (cell)  files 

MOSS  (vector)  files  to  MAPS  (cell)  files 

ADS  Rev.6  files  to  ADS  Rev.7 

MAPS  raster  data  to  ASCH  file. 

Point,  line  and  polygon  data  to 
MOSS  IMPORT/EXPORT  format. 


ADS  menu 
ADS  menu 
ADS  EDIT  option 

IS>ADS>CNVADS.RUN 
MOSS  command 
MOSS  UTILITY 
ADS  menu 

IS>MOSS>DLGS20.RUN 
MOSSUnLTTY 

MOSS  command 
ADS  menu 
MOSS  command 
MAPS  command 

MOSS  command 
CMDNCOCPL 
MOSS  command 
IS>MOSS>MOSS2LGS.RUN 
ADS  menu 

ADS  menu 
ADS  (stand  alone) 

CMDNCOCPL 

MOSS  command 
MAPS  commanc 
ADS  menu 
MOSS  command 
MOSS  UTILITY 
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PRIMOS  Commands 


The  PRIMOS  commands  most  frequently  used  during  GIS  sessions  are  listed 
below.  For  more  detailed  information,  see  the  PRIMOS  Commands  Reference 
Guide  or  use  the  PRIMOS  HELP  command. 


AB-AC 

ABBREV 

ATTACH 

CHANGE_PASSWORD 

CLOSE  ALL 

CNAME 

COMOUTTUT 

CONTROL P 

COPY 

DELETE 

DGWN 

EDIT_ACCESS 

EMACS 

HELP 

ICE 

LD 

UST_ACCESS 

LOGIN 

LOGOUT 

SUST 

SPOOL 

STAT^'*^ 

SYSTEM 

UP 

WHERE 


creates  abbreviations  in  the  text  editor 

defines  user's  abbreviations 

changes  user's  current  directory 

changes  a  user  login  password 

closes  all  files  that  may  be  opened 

changes  a  file  name 

opens  "como"  file  and  stores  log  of  session 

terminates  current  processing 

copies  file  or  directory 

deletes  a  file  or  directory 

changes  to  directory  below  current  directory 

edits  existing  access  control  hst* 

accesses  the  text  editor 

provides  information  about  PRIMOS  commands 

initiates  logout/login  automatically 

lists  contents  of  directory 

lists  access  of  a  file  or  directory* 

initiates  PRIMOS  session 

terminates  PRIMOS  session 

displays  file  contents  on  screen 

sends  files  to  line  printer 

displays  system  information  regarding  disks,  users,  and  depict 

executes  PRIMOS  command  with  MOSS 

changes  to  directory  above  current  directory 

gives  pathname  of  current  directory 


*  Access  codes  are:  Protect;  Delete;  Add;  List;  Use;  Read;  Write,-  ALL;  and  NONE 
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Disk  Management 

Disk  Space 


The  amount  of  disk  space  varies  with  the  type  of  data  used.  For  example,  one 
public  land  survey  map  required  55,000  bytes  (27  records)  in  vector  form  and 
78,000  bytes  (38  records)  in  raster  (cell)  form.  (One  record  equals  2,048  bytes 
or  16,384  bits.) 

Disk  space  also  varies  with  thematic  types.  A  few  data  themes  and  their 
estimated  disk  space  requirements  are  as  follows: 

THEME  Megabytes  per  7-1/2  minute  quadrangle 


LANDNET 

0.07 

GEOLOGIC 

0.557 

SOIL 

0.548 

SURFACE  WATER 

0.021 

ELEVATION 

4.35 

When  you  are  determining  space  requirements  for  MFDs  and  UFDs,  the 
following  PRIMOS  commands  can  be  useful: 

AVAIL      Displays  disk  usage  statistics  (total  number  of  records  and 
number  of  available  records  in  a  directory) 

SIZE  Displays  the  size  of  a  file  or  the  number  of  entries  in 

a  directory 

ID  Displays  the  number  of  records  and  quota  for  a  directory, 

along  with  the  directory  contents 

*  1  record  =  2048  bytes  or  16,384  bits 

File  Cleanup  and  DEiEnoNS 

In  PRIMOS,  use  DELETE  to  delete  a  file  -  r  directory. 

In  MOSS/MAPS,  use  DELETE  to  delete  maps  from  yotir  working  directory. 
Maps  in  the  master  directory  can  only  be  deleted  if  you  have  access  to  do  so. 
The  DELETE  command  will  erase  the  map  file  from  yoiu:  directory  and 
remove  the  name  of  the  map  from  your  polygon.dt  file.  If  you  want  to ' 
move  the  map  name  from  the  polygon.dt  file  but  not  erase  the  map  fil  . 
use  the  command  UTILITY-DATABTEST  option. 

Archiving 

In  order  to  free  space  on  the  system  disk  for  other  processing,  contact  your 
systems  operator  to  archive  inactive  data  and  map  files  onto  a  magnetic  tape. 
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The  Text  Editor  (EMACS) 


You  can  use  tlie  EMACS  text-editing  program  to  create  and  edit  text  files. 
For  example,  you  can  use  it  to  create  CPL  (Command  Processing  Language) 
files  to  be  used  for  batch  jobs  and  data  files  to  be  used  as  multiple  attribute 
files  in  MOSS. 

Before  using  EMACS,  you  must  use  the  reset  command  to  set  the  termuial 
CODE  to  ANSI.  To  invoke  EMACS,  you  must  use  the  command 

EMACS  -TTP  VTIOO  filename  -NOXOPP 

Usmg  the  PRIMOS  command  "AB-AC",  you  can  create  a  two-letter  abbrevia- 
tion that  invokes  the  EMACS  editor.  It  is  much  easier  to  type  a  two-letter 
abbreviation  than  to  type  the  long  command  every  time  you  want  to  use 
EMACS.  The  following  example  creates  an  abbreviation  "EM": 

AB  -AC  EM  EMACS  -TTP  VTIOO  %1%  -NOXOFF 

A  brief  Usting  of  useful  EMACS  commands  are  listed  in  Table  1.2. 
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TABLE  1.2.  Useful  EMACS  commands 


Viewing  Screens: 

CTRLV 

ESCV 

CTRLL 

Cursor  Movement: 

CTRLF 

CTRLB 

ESCF 

ESCB 

ESCE 

ESC  A 

CTRL] 

CTRL! 

CTRLN 

CTRLZ 

CTRLA 

CTRLE 

ESC< 

ESC> 

Saving  and  Creating  Files: 

CTRLXCTRLF 

CTRLXCTRLL 

CTRLXCTRLS 

CTRLXCTRLW 

CTRLXCTRLC 

Removing  and  Restoring  Text: 

(bksp) 

CTRLD 

ESCD 

ESC  (bksp) 

CTRLXCTRLK 

CTRLK 

CTRLY 

ESCY 

CTRLW 

Marking  Text  (The  cursor  and  mark  specify  the  ends  of  a  region.) 

CTRL  @  Set  a  mark  at  cursor 

CTRL  X  H  Mark  entire  buffer 

CTRL  X  CTRL  X  Exchange  the  positions  of  the  mark  and  cursor 

Searching  and  Replacing: 

CTRLS 
CTRLR 
ESC  X  Replace 

ESC  % 


Move  forward  one  screen 

Move  backward  one  screen 

'Refresh'  the  current  screen  and  center  screen  where  cursor  is 

Move  forward  a  character 
Move  backward  a  character 
Move  forward  one  word 
Move  backward  one  word 
Move  forward  one  sentence 
Move  backward  one  sentence 
Move  forward  one  paragraph 
Move  backward  one  paragraph 
Move  down  one  line 
Move  up  one  line 
Move  to  beginning  of  line 
Move  to  end  of  line 
Go  to  beginning  of  file 
Go  to  end  of  file 

Find  or  create  file  you  specify 

Insert  a  file 

Save  changes  in  current  file  (some  systems  use  CTRL  X  ESC  S) 

Write  to  a  file 

Exit  EMACS  and  return  to  PRIMOS 

Delete  the  character  immediately  before  cursor 

Delete  the  character  immediately  after  cursor 

Delete  next  word 

Delete  previous  word  ~ 

Delete  text  from  start  of  line  to  cursor 

Delete  text  from  cursor  to  end  of  line  and  put  in  kill  ring 

Yank  back  text  in  kill  ring  and  insert  it  in  file  at  cursor 

Yank  prior  text  from  kill  ring  and  replace  most  recently  yanked  item 

Delete  region  bounded  by  the  mark  and  cursor  and  copy  to  kill  ring 


Help  Commands: 

CTRL_ 

C 
A 
D 
L 
? 


Search  forward  (some  systems  use  ESC  S) 
Search  reverse 

Globally  replace  one  string  with  another 
Query  before  replacing 

Invoke  online  help 

Explain  the  command  typed  in  by  user 

List  commands  related  to  subject  typed  in  by  user 

Invoke  the  "Describe"  function 

List  last  20  characters  typed 

Summarize  all  of  above  options  at  top  of  =creen 
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TABLE  1 .2.  Useful  EMACS  commands  (continued) 


Capitalize  word 

Make  work  all  upper  case  word 

Make  work  all  lower  case  word 


Case  Conversions 

ESCC 

ESCU 

ESCL 

Macros: 

CTRLXI 

CTRLX) 

CTRLXE 

ESC#CTRLC 

Other  Useful  Commands 

CTRLG 

CTRLO 

ESC# 

CTRLC 

CTRLT 

ESCT 

Begin  macro  definition 

End  macro  definition 

Execute  most  recently  defined  keyboard  macro 

Repeat  the  last  command  #  times  (use  CTRL  X  E  first  for  a  macro) 


Cancel  partially  typed  command 
Insert  new  Une  text 
Repeat  next  command  #  times 
Repeat  previous  action  or  command 
Transpose  characters 
Transpose  words 
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Non-Interactive  Processing 

Non-interactive  processing  allows  you  to  run  a  series  of  commands  within  a 
CPL  file.  The  advantages  of  non-interactive  processing  are  that  you  can: 

*  Continue  other  tasks  without  being  tied  to  the  display  terminal 

*  Create  long,  complex  strings  of  commands  or  processes  that  you  can  use 
again 

*  Use  the  CPL  file  to  log  processes  for  a  project 

*  Easily  alter  command  parameters  or  variables 
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Create  a  CPL  File 

Refer  to  the  following  example  of  a  CPL  file  named  TEST.CPL.  Comment 
lines  begin  with  /*.  The  nimiber  of  CPL  lines  for  a  MOSS  command  must 
be  the  same  as  if  you  were  using  the  command  interactively. 

/*   This   is   a  MOSS  processing  CPL  named  TEST.CPL 


&riATA  MOSS 

MINTRAIN 

TEAMl 

SYSTEM    COMO    TEST.ODT 

OPEN    SJCOAL 

FREE    ALL 

SELECT    LANDLINES    ALL 

WINDOW    1 

SELECT  COAL.PTS  ALL 

SELECT  S J. LEASE  ALL 

ACTIVE 

PLOT  2  3 

BYE 

&END 

CCadO  -END 

&RETDRN 


/*  invokes  MOSS  iintil  &END 

/*  MOSS  project  directory 

/*  MOSS  user  naaa 

/*   opens  como  file  "TEST. OUT" 

/*  opens  desired  MOSS  project 

/*  frees  existing  active  table 

/*  select  map 

/*  window  to  this  map 

/*  select  second  map 

/*  select  third  map 

/*  list  active  table 

/*  plot  maps  2  and  3 

/*  exit  MOSS  session 

/*  end  of  MOSS  data 

/*  close  como  file 

/*  end  of  CPL  file 
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Submit  a  CPL  file  as  a  batch  job 

To  submit  a  CPL  file  as  a  batch  job,  simply  enter  the  PRIMOS  command: 
JOB   filename 

After  you  have  submitted  the  job  to  the  batch  queue,  the  batch  Job 
Identification  Number  will  be  displayed. 

Controlling  a  batch  job 

These  PRIMOS  commands  can  be  used  to  examine  or  modify  the  status 

of  an  existing  batch  job: 

JOB  #  -ABORT  Terminates  a  running  job 

JOB  #  -CANCEL         Cancels  a  waiting  job 
JOB  #  -CHANGE         Changes  entry  characteristics 
JOB  #  -STATUS         Shows  a  job's  progress 

Phantoms 

Another  way  to  execute  a  list  of  commands  in  a  non-interactive  processing 
mode  is  to  use  phantoms.  Unlike  PRIMOS  batch  jobs,  phantoms  allow  CPL 
files  to  run  immediately  without  being  submitted  to  a  queue.  To  start  a 
phantom,  enter: 

PHANTOM      filename. CPL 

The  system  will  then  display  a  phantom  user  number. 

If  a  phantom  is  used  to  run  MOSS/MAPS,  you  cannot  enter  MOSS/MAPS 
unless  you  use  a  different  directory.  However,  you  can  do  work  at  the  PRI- 
MOS level  in  the  same  directory  while  nmning  a  phantom. 

If  you  want  to  cancel  a  phantom,  use  the  STATUS  ME  command  to  list  the 
phantom  user  nimiber.  Then  enter: 

LOGOUT      -phantom  user  number 
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GIS  Project  Management 


The  following  project  management  checklist  is  a  general  procedure  common 
to  many  GIS  projects.  Some  projects  may  require  more  detailed  steps  in 
certain  areas.  For  example,  a  project  using  remote  sensing  data  may  need 
more  specific  plaiming  in  data  acquisition  and  classification. 

Project  Management  •-  A  Checklist 

I.  Define  project 

A.  Objectives  and  altematives 

B.  Extent  of  study  area 

-  representative  area 

-  manageable  size 

-  available  data 

C.  Plan  of  attack 

-  time,  money,  personnel  constraints 

-  schedule;  deadlines 

-  milestones  and  interim  products 

-  responsibilities  and  outside  contracting 

-  final  products  and  their  uses 

D.  Capabilities  of  system 

E.  Management  priority  areas 
n.  Define  Data  Requirements 

A.  Themes  (one  theme  per  map) 

-  coal  data  points 
-streams 
-roads 

-  landlines 

-  ownership 

-  lease  boundaries 

B.  Required  data  types 

-point 
-line 

-  polygon 

-  DEM  (Digital  Elevation  Models) 

-  DLG  (Digital  Line  Graphs) 

-  cell  maps  (dichotomous,  discrete,  or  continuous) 

C.  Data  considerations 

-  accuracy 

-  availabihty 
-vintage 

-  resolution 

-  scale  and  projection  compatibihties 
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D.  Data  management 

-  map  naming  conventions 

-  data  backups 

-  archiving 

-  tape  library 

-  security 
in.  Acquire  data 

A.  Previously  digitized  data 

B.  Other  agencies  as  sources 

-  Bureau  of  Mines:  Mining  Industry  Location  System  (MILS)  data 

-  Geological  Survey. 

Mineral  Resources  Data  System  (IVIRDS)  data 

DEM 

DLG 

-  Defense  Mapping  Agency:  DEM 

C.  Data  acquisition 

-  field  work 

-  air  photo 

-  satellite  digital  images 

IV.  Determine  system  use 

A.  Storage  capacity 

B.  Processing  time 

C.  Hardware  maintenance  and  upgrades 

D.  Plaimed  and  unplanned  downtime 

E.  Software  enhancements 

F.  Scheduling  access 

G.  Project  directory  management 

V.  Prepare  data 

A.  Desired  interim  and  final  products  according  to  software  capabilities 

B.  Accuracy  and  missing  data  checks 

C.  Format  requirements;  reformatting  if  necessary 

D.  Adjacent  maps  for  later  edge-fitting  requirements 

E.  QuaUty  control 

F.  Time  and  money  comparison 

G.  Compilation  sheets 

-  stable  media 

-  common  errors: 

edge  matching 
inconsistent  labeling 
missing  data 
extraneous  data 
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VI.  Enter  data 

A.  Loading  previously  digitized  and  acquired  data 

B.  Prioritizing  remaining  data  entry  according  to 

-  most  urgent  need 

-  data  complexity 

-  data  availability 

C.  Digitizing  steps 

(approximately  3  hours  per  7-1/2  minute  quadrangle  map) 

-  register  map  to  tablet 

-  digitize  in  ADS 

-  plot  resultant  map  and  compare  to  original,  i.e.,  verify 

-  edit  and/or  digitize  to  correct  errors  and  add  data  plot 
and  verify 

D.  Documenting  maps  and  themes  entered  via  naming  conventions 

E.  Managing  newly  entered  data 

-  map  names 

-  project  files 

-  archiving 
-security 

Vn.  Conduct  data  transfers,  reformats,  and  system  interfaces  (See  Figure  1.1) 

A.  ADS  to  MOSS  (ADS2MOSS) 

B.  MOSS  to  ADS  (MOSS2ADS) 

C.  DEM  to  MAPS  (IMPORT  command  in  MAPS) 

D.  DLG  to  MAPS  (DLGS20,  DLG02ADS  or  MOSS  Utility 
command-DLG3  option) 

E.  MILS  data  to  MOSS  import  format 

F.  MRDS  data  to  MOSS  import  format 
Vin.  Analyze  data 

A.  Project  objectives  and  choice  of  analyses 

B.  Construction  of  models 

C.  Interim  products 
IX.  Output  products 

A.  Limitations  of  output  capabilities 

-  resolution 

-  accuracy 

-  aesthetics 

-  appropriateness 

B.  Tables,  graphs,  and  text  output 

-  text  editor  (EMACS) 

-  MOSS/MAPS  text  output 
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C.  Maps 

-  COS  (Cartographic  Output  System)  maps 

-  cartographic  design  of  output  map 

scaleshade  and  color 
theme(s)legend 
complexity  lettering 
symbolization 

-  MOSS  PENPLOT  command  (for  vector  maps) 

-  MOSS  GENPLOT  command  (for  cell  and  vector  maps) 

-  ADS  production  plots 

-  GISPLOT,  which  reformats  files  from  COS,  PENPLOT, 
and  GENPLOT,  and  creates  a  plot  file  for  specified  plotter 
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PRIMOS  Review  Exercise  1 

The  purpose  of  this  exercise  is  to  review  selected  PRIMOS  commands. 

1 .  LOGIN  to  the  PRIME  using  your  user  ID  and  password. 

2.  Go  down  (DOWN)  to  the  MOSS  directory.  Find  your  location  in  the 
system  (WHERE  command).  List  (LD)  the  contents  of  the  directory. 
Go  back  up  to  your  top  level  directory  (UP  or  ORIGIN). 

3.  Use  the  HELP  command.  Respond  with  NO  when  you've  seen  enough 
of  the  Usting  of  cormnands.  Use  the  HELP  command  to  obtain  informa- 
tion on  the  ABBREV  command 

4.  List  (LD)  the  files  in  your  directory  to  see  if  you  have  an  ABBREV  file 
already  created.  List  the  abbreviations  in  the  ABBREV  file. 

If  there  is  no  ABBREV  file,  create  one  using  the  ABBREV  -CR  command. 
Activate  the  new  ABBREV  file  by  using  EIVIACS  to  add  the  following 
Une  to  your  login.cpl: 

ab    filename . abbrev 

Do  an  ICE  (Initialize_Activate_Enviroimient)  to  activate  the  ABBREV  file. 

5.  Add  the  following  abbreviations  to  your  ABBREV  file: 

W  where 

Q  close  -alljrls  -alljcomo  -end 

DEL  delete 

SL  slist 

6.  Set  up  the  abbreviations  to  start  MOSS  with  a  single  keystroke. 

AB    -AC   M  MOSS   HINTRAIN  usemame 

7.  LAC  (List_Access!  to  display  the  ACL's  (Access_ControLList)  of  a  file 
in  your  directory. 

8.  Use  SIZE  command  to  list  file  and  directory  sizes. 

9.  Use  STATUS  US  and  STATUS  ME  to  display  user  information. 

10.  Create  and  run  a  CPL  file  that  starts  a  como  file,  selects  and  pi  .3 
several  maps,  lists  the  active  table,  frees  that  active  table  and  then 
closes  the  como  file. 
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Chapter  Two: 

GIS  Applications 
for  Coal  Resource 
Analyses 

In  this  chapter  you  will 

■  Generate  a  boundary  map 
for  the  study  area 

u  Generate  a  coal  thickness  grid 

■  Evaluate  and  customize  a  grid 

u  Calculate  Reserve  Base  coal  tonnage 
m  Generate  a  stripping  mining-ratio  grid 

■  Generate  a  mining-ratio  grid  using 
multiple  Recovery  Factors 

■  Create  a  coal  potential  map 

m  Generate  a  suiface-minable  coal  grid 

B  Prepare  a  DEM  for  map  analyses 

m  Compare  coUar  elevations  to  the  r/M 

m   Use  DEM  to  generate  an 
overhr-'den  grid 
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GIS  for  Coal  Resource  Analysis  —  introduction 

A  Geographical  Information  System  (GIS)  is  well-suited  for  coal  resource 
assessment  activities,  such  as  generating  isopach,  structure  and  overburden 
maps,  calculating  reserves  and  ranking  mining  potential  of  coal  deposits.  All 
the  information  in  a  GIS  data  base  is  identified  geographically,  that  is,  by  X 
and  Y  coordinates.  These  coordinates  can  be  in  degrees  (latitude/longitude), 
meters  (UTMs)  or  feet  (State  Plane),  but  must  be  consistent  within  a  study 
area.  Attached  to  this  geographically  located  information  are  data  elements, 
or  attributes.  These  attributes  can  include  Z- values  (e.g.,  ground  surface 
elevation,  depth  to  a  coal  seam,  thickness,  analytical  data)  and  general  infor- 
mation (e.g.,  drill  hole  identifier,  date  of  sampUng). 

Once  the  data  are  in  the  GIS  environment,  the  Z- values  can  be  transformed 
into  gridded  numerical  surfaces.  These  grids  (cell  maps)  can  then  be  added, 
subtracted,  multipHed,  divided,  or  algebraically  manipulated.  For  example, 
you  can  generate  an  overburden  grid  by  subtracting  an  elevation  grid  for  the 
top  of  a  coal  seam  from  the  surface  elevation  grid.  Keep  in  mind  that  these 
spatial  analyses  can  only  be  done  with  grids  that  have  the  same  number  of 
rows  and  columns  and  the  same  cell  size;  the  accuracy  of  the  GIS  cannot  be 
any  greater  than  the  data  you  put  into  the  data  base.  New  data  available  can 
easily  be  added  to  the  data  base  to  generate  new  grids. 

Contour  values  for  thickness,  bed  structure,  overburden  or  stripping  mining- 
ratios  may  extend  beyond  the  area  of  interest.  This  may  include  botmdaries 
(lease,  mined  areas,  bum  areas,  channels)  or  other  coal  limitations  (chemical, 
minimum  thickness).  In  MOSS,  the  contour  hnes  can  be  truncated  at  the 
polygon  boundaries  using  the  LPOVER  command.  Lines  of  limitation  (e.g., 
the  zero  thickness  contour  line)  must  first  be  converted  to  a  polygon  map 
using  the  M0SS2ADS  and  ADS  commands.  In  MAPS,  the  cell  values  can 
be  limited  to  the  area  of  interest  by  using  the  RENU2VIBER  and  MATH 
commands. 

You  can  generate  customized  maps  tl:    .  will  show  the  maps  as  overlays. 
This  can  also  include  other  maps  such  is  transportation  networks,  streams, 
political  boundaries,  leases,  townships  and  land/mineral  ownership.  The 
ways  in  which  they  can  be  displayed  are  as  follows: 

a.  Use  the  display  commands  such  as  PLOT,  FLOOD,  and  SHADE  to  view 
the  maps/grids  on  the  terminal  screen. 

b.  For  vector  maps  only: 

(1)  Use  the  PENPLOT  command  in  MOSS  to  create  2  metafile,  which  can 
then  be  plotted  on  the  screen  or  a  plotter  using  GI£?LOT  in  PRIMOS. 
These  maps  can  be  plotted  to  the  map  scale  of  your  choice,  with  the 
size  of  the  plotter  paper  being  a  limitingfactor. 
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(2)  Use  the  DISPLAY  command  in  MOSS  to  build  a  display  file 
(filenameDP),  which  can  be  viewed  on  the  terminal  screen  using 
SHOW  command  in  PRIMOS. 

c.  For  cell  maps  only: 

(1)  Use  the  GENPLOT  command  in  MOSS  to  create  a  metafile,  which  can 
be  plotted  on  the  screen  or  a  plotter  using  GISPLOT.  The  maps  can  be 
plotted  to  the  map  scale  of  your  choice,  with  the  size  of  the  plotter 
paper  being  a  limiting  factor. 

(2)  Use  the  DISPLAY  command  in  MAPS  to  create  a  metafile  or  a  display 
file.  The  prefix  "M."  is  automatically  added  to  the  metafile;  the  suffix 
".DP"  is  automatically  added  to  the  display  file.  The  metafile  can  be 
plotted  on  the  screen  or  a  plotter  using  GISPLOT  in  PRIMOS.  The 
display  file  can  be  seen  on  the  terminal  screen  using  the  VIEW  com- 
mand in  MAPS. 

d.  Use  the  3D  command  to  drape  cell  maps  over  an  elevation  model. 

You  can  add  text  using  the  TEXT  or  ANNOTATION  commands  and  label 
contours  using  the  AUTOLABEL  or  LEGEND-CONTOUR  commands.  The 
Cartographic  Output  System  (COS)  allows  you  to  produce  high  quality  maps 
by  using  Bureau  mapping  standards  for  fine  styles,  symbology,  shading 
pattems,  and  textual  information. 

Flow  Chart  2.1  shows  the  CIS  process  for  a  coal  resource  assessment  of  a 
study  area.  It  was  necessary  to  use  two  study  areas  to  take  you  through  most 
of  the  map  analyses  shovm  in  the  flow  chart.  The  first  study  area  has  no 
Digital  Elevation  Model  (DEM).  Thus  map  analyses  with  the  DEM  are 
provided  by  a  second  study  area. 

Using  a  CIS  may  lead  to  new  approaches  to  coal  resource  management, 
because  it  enables  the  geologist  to  quickly  evaluate  a  combination  of  maps 
and  study  a  problem  in  ways  that  would  not  be  possible  if  these  maps  had 
to  be  compiled  in  a  conventional  manner.  Geographical  relationships  of 
multiple  themes  that  can  facilitate  the  decision-making  process  can  easily 
be  shown. 
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Flow  Chart  2.1 .    Processes  used  for  map  analyses 

Part  A:  Generate  the  thickness  grid 
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Flow  Chart  2.1.  (continued) 

Part  B:  Calculate  Reserve  Base  tonnage 
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Flow  Chart  2.1 .  (continued) 

PartC: 

Use  DEM  to  Generate  Overburden  and  Stripping  Mining-Ratios  Grids 
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Flow  Chart  2.1.  (continued) 

PartD: 

Generate  Surface-Minable  Coal  tonnage 


StriniingMtalBg-ilatioGrid 
(siiigia  m  wiltipic  itacvniy  FmIms) 

CMtiMMd 


Generate  Surface-Mimble  Coal  (SMC)  grM 


•  Cron  ration  grid  wMh 

■iiiafale<eal  grid.  assigB  value  1 

to  celb  within  UbHs 


Generate  ceal-llBiUtioa  gridi 


•  bid  analyticai  vaiuea 
(aifc,  sulfur,  etc.) 

|_     _     _    «Ml_     —     — 

•  INgHiza  and  rasteria 
etkercaaiiieiitatioas 

(bun  or  ■Ined  areas,  duMwIs) 


Generate  ■taable<eal  grid 


•  Cress  grids,  assign  values  1  to 
eeilswitfiiniiaiU 


Generate  SIWC  grid  contained  witliin  lease 


•  Cross  SMC  grid 
witfe  rasterized  lease 


I 


•  Label  SMC-lease  grid; 
Maintain  ceil  values; 
enter  subject  naaos 


I 


•  Total  thickness  grid  by 
SMC^ease  grid  for  tonnage  rapoit 


not  eannd  in  ttiit  tnkiing  ■amrai 


CHAPTER  TWO:  QIS  Apalications  for  Goal  Rssourc*  AmaSyscs 


31 


Given  Data  for  the  SJCO AL  study  area 

The  SJCOAL  study  area  is  located  in  New  Mexico.  Six  townships  are  in- 
volved in  the  SJCOAL  study  area:  T23-25N  and  Rl  1-12W.  There  are  45  coal 
drill  holes  in  the  western  part  of  the  study  area. 

Master  maps  in  the  MOSS  project  SJCOAL: 

IVIaster  Map  Name        Map  Type  Map  Theme 


COALJTS 

1  (potQt) 

Coal  data  for  45  drill  holes 

COAL.SAMP 

1  (point) 

Coal  sample  data  (ash,  sulfur,  etc.) 

LANDT,TNES 

3  (polygon) 

Township,  range,  section 

SJ.LEASE 

3  (polygon) 

Coal  lease  boundary 

ROADS 

2  (line) 

Roads 

STREAMS 

2(Une) 

Streams 

WSA 

3  (polygon) 

Wilderness  Study  Areas 
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COAL  EXERCISE  1 

Familianze  yourself  with  the  master  maps. 

In  this  exericise  you  will  list,  select  and  plot  the  master  maps  in  order  to 
become  familiar  with  the  data.  You  will  learn  how  to  do  a  Boolean  search  to 
retrieve  data  and  how  to  calculate  the  elevation  at  the  top  of  the  coal  bed 

a.  Enter  MOSS  and  the  MOSSDATA  directory  called  MINTRAIN. 
OPEN  the  project  SJCOAL. 

b.  LIST  the  maps  ia  the  master  directory  for  this  study  area. 

c.  Use  UmiTY-BROWZE  to  list  the  master  maps  and  their  themes. 

The  output  file  will  automatically  be  named  SJCOAL.BL.  Use  the  SIIST 
command  on  the  output  file  (in  MOSS,  use  SYSTEM  SUST  SJCOAL.BL). 

d.  There  is  a  map  named  COAL.PTS  which  contains  45  drill  sites  and 
associated  attributes. 

(1)  DESCRIBE  this  map  to  become  familiar  with  the  multiple  attributes, 
such  as  coal  thickness,  overburden,  interburden,  etc. 

(2)  SELECT  and  WINDOW  to  COAL.PTS.  PLOT  and  QUERY  the  mapto 
obtain  information  about  the  attributes. 

(3)  Create  a  REPORT  with  a  tabular  listing  of  several  attributes  associated 
with  the  drill  sites. 

e.  To  find  areas  that  have  a  low  mining  ratio,  use  the  BSEARCH  command 
to  do  a  Boolean  search-and-select  from  the  COAL.PTS  map  to  find  those 
drill  holes  in  which  the  overburden  is  less  than  30  feet  and  total  coal 
thickness  is  greater  than  15  feet.  PLOT  the  selected  drill  holes  and  use  the 
LEGEND  command  to  label  these  points  by  their  drill  hole  ED. 

f.  SELECT  and  WINDOW  to  the  township  and  range  map  (LANDLINES). 
SELECT  the  rest  of  the  master  maps,  excluding  the  subject  "null"  from 
the  WSA  map. 

FLOOD  the  lease  map  (SJ.LEASE)  in  orange 

SHADE  the  Wilderness  Study  Areas 

PLOT  the  landlines  in  white 

PLOT  the  streams  in  dark  blue 

PLOT  roads  in  red 

PLOT  the  data  points 


?H&i>TER  T^O:  ©IS  Ajiplicatlsiis  ior  C««l  R«seurc«  Analyses       33 


COAL  EXERCISE  1  (cont'd) 


g.  The  COAL.SAMP  map  in  the  data  base  contains  infonnation  on  moisture, 
ash,  carbon,  BTU's,  and  other  sample  data. 

(1)  DESCRIBE  the  attributes  for  this  map. 

(2)  Use  the  BSEARCH  command  to  select  the  sampled  coal  holes  which 
have  a  sulfur  content  greater  than  0.5  percent. 

(3)  PLOT  the  selected  coal  holes.  Label  (LEGEND)  the  point  locations  by 
sulfur  content. 

h.  Create  a  new  attribute: 

(1)  SAVE  the  COAL.PTS  map  and  name  it  MYCOAL.PTS  since  you  do 
not  have  write  permission  to  modify  a  master  map. 

(2)  SELECT  MYCOAL.PTS 

(3)  Use  the  COMPUTE  command  to  calculate  the  elevation  at  the  top  of 
the  coal  bed  (elevation  minus  total  overburden). 

*  Name  the  new  attribute  "coai.top" 

(4)  Generate  a  report  that  lists  the  elevation,  total  overburden  and  the  new 
attribute  "coai.top". 
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Flow  Chart  2.2  Generate  a  boundary  map  for  the  study  area 


SELECT,  WINDOW  and  PLOT 
point  and  lease  maps 


T 


ZOOM,  POINT  option  to 
ciiange  window  size 


I 


ERASE,  rePLOT  maps 


GENERATE 
boundary  map 


I 


SELECT  and  WINDOW  to 
boundary  map 
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COAL  EXERCISE  2.  Generate  a  boundary  map  for  the  study  area 

Cieating  a  boundary  map  for  the  study  area: 

a  maxmiizes  screen  view  of  the  data  points 

n  Bmits  interpolation  to  area  that  is  controlled  by 
data  points 


a.  GENERATE  a  boundary  map  that  includes  all  the  points  in  the 
MYCOAL.PTSmap: 

(1)  SELECT  the  points  map  (MYCOAL.PTS)  and  the  lease  map  (SJ.LEASE). 
WINDOW  to  both  maps  and  PLOT  them. 

(2)  Because  some  points  are  at  the  edge  of  the  map,  use  the 
following  procedure: 

(a)  Use  ZOOM  POINT  and  a  magnification  factor  of  0.8  to  change 
the  window.  Place  the  crosshairs  in  the  center  of  the  map  on 
the  screen. 

(b)  ERASE  the  screen  and  rePLOT  the  maps.  All  points  should  now 
be  within  the  edges  of  the  map. 

(3)  Using  the  zoomed  window,  GENERATE  a  boundary  that  includes  all 
the  points  and  the  lease  boundary. 

*  Name  the  new  map  SJ.BORDER. 


b.  SELECT  and  WINDOW  the  SJ.BORDER  map.  PLOT  the  points  map 
(MYCOAL.PTS)  and  the  lease  map  (SJ.LEASE).  All  the  points  should 
now  be  within  the  border. 
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The  ABCs  of  Gridding 

The  Gridding  Process 


The  gridding  process  generates  a  raster  map  over  the  study  area  and,  using 
one  of  the  gridding  algorithms,  fills  in  the  probable  value  for  each  cell.  The 
process  looks  at  points  of  observed  values  and  then  calculates  values  for  the 
void  cells  where  measvirements  have  not  been  taken.  The  resultant  grid  can 
be  used  to  generate  a  contour  map  (e.g.,  isopach  map)  or  to  conduct  further 
analyses  (e.g.,  subtracting  structure  grid  from  the  elevation  grid  to  produce 
an  overburden  grid). 

The  GRID  command  in  MOSS  provides  five  different  gridding  algorithms  for 
determining  the  values  for  each  cell.  Three  of  these  algorithms  are  illustrated 
in  Figure  2.1.  Also  see  Table  2. 1 .  for  a  comparison  of  the  advantages  and 
disadvantages  of  the  various  gridding  algorithms. 


Four-Point  Quadrant  Whghted  Average 

The  FOUR-POINT  method  uses  each  cell  as  the  center  of  a  quadrant.  The 
algorithm  determines  the  closest  data  point  to  the  center  of  the  ceil  in  each 
of  the  four  quarters  and  calculates  their  weighted  average  to  determine  the 
value  (z)  for  that  cell.  The  greater  the  distance  of  the  data  point  from  the  cell, 
the  less  influence  the  data  point  has  on  determining  the  value  of  the  cell.  It 
then  moves  on  to  the  next  cell  until  all  cells  within  the  grid  have  a  value. 
The  value  (z)  of  the  cell  is  determined  by  this  formula: 

2=   zl/dl  -I- z2/d2  +  z3/d3  +  z4/d4 
1/dl  +  l/d2  +  l/d3  +  l/d4 

Eight-Point  Quadrant  Weighted  Average 

For  each  cell,  the  EIGHT-POINT  method  searches  for  the  nearest  poinii  that 
are  within  the  specified  moving-window  size.  It  must  find  at  least  3  points 
within  the  window  and  limits  itself  to  a  maximimi  of  8  nearest  points.  Any 
cell  with  less  than  3  nearest  points  within  the  window  is  assigned  a  value  of 
10  less  than  the  minimum  value  of  the  point  map.  A  weighted  average  is 
computed  based  on  the  distances  from  the  ceil  to  each  point.  The  value  iz! 
of  a  cell  is  determined  by  the  same  formula  as  the  FOUR-POINT  m*   .  •':: 
except  eight  point- values  are  used  rather  than  four. 
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Simple  or  Universal  Krigevg 

Kriging  is  a  moving-average  method  that  uses  the  theory  of  regional  variables 
(e.g.,  thickness  of  a  coal  bed).  To  determine  input  values  for  kriging,  the 
VARIOGRAM  command  must  first  be  used  to  generate  a  semi-variogram. 
You  should  not  attempt  to  run  the  kriging  option  without  some  understand- 
ing of  the  theory  of  geostatistics. 

QuiNnc  Spline  iNTERPRETAnoN 

The  spUne  method  divides  the  area  into  a  network  of  triangles  using  the 
original  data  points  as  vertices.  It  then  uses  the  nearest  3-10  data  points 
(determined  by  user  input)  to  estimate  slope  and  gradient  at  each  data  point. 
The  slope  and  gradient  is  used  to  calculate  a  fifth  degree  (quintic)  interpolat- 
ing polynomial  over  each  triangle.  This  causes  adjacent  triangles  to  match 
smoothly  (to  spline).  Finally,  the  program  determines  in  which  triangle  the 
center  of  the  cell  Ues  and  the  cell  value  is  calculated  based  on  the  polynomial 
over  this  triangle. 

Nearest  Nhghbor  (Polygon) 

The  nearest  neighbor  method  determines  the  closest  data  point  to  the  target 
cell  and  assigns  the  target  cell  that  value.  If  more  than  one  data  point  are 
eqtially  closes  to  the  target  cell,  the  cell  is  assigned  the  average  of  the  values. 
The  resulting  ceil  map  appears  to  be  made  up  of  polygons. 

CTOG  (Contour  to  Grid) 

CTOG  is  a  MOSS  command  that  allows  you  to  create  a  grid  (cell  map)  from 
an  existing  contour  map  (line  map.  Type  2).  You  must  first  create  an  export 
file  firom  the  contour  map  in  MOSS  using  the  EXPORT  command.  Secondly, 
use  the  UTILITY-CTOG  option  to  reformat  the  export  file  into  a  MAPS 
import  file.  Thirdly,  use  the  IMPORT  command  in  MAPS  to  generate  a 
grid  with  the  import  file,  using  the  "DEM  7.5"  format  option. 
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Factors  to  remember  when  gridding 

You  must  set  the  window. 

The  Minimum  Boimding  Rectangle  (MBR)  of  the  resultant  grid  may  be 
slightly  different  than  the  original  window. 

If  the  grid  is  used  in  the  volumetric  calculations,  in  most  cases  the  acreage 
will  be  different  than  the  vector  map  acreage.  Thus  use  the  average  thickness 
obtained  through  this  method  and  then  manually  calculate  the  tonnage  using 
the  acreage  from  the  vector  map. 

You  cannot  grid  with  dupUcate  points  in  the  Quintic  Spline  method.  Use  the 
Variogram  command  to  determine  whether  you  have  duplicate  points  in  your 
data  set.  If  so,  remove  the  dupUcates  before  gridding. 

If  a  very  small  cell  size  is  chosen  and  number  of  rows  and  columns  are  large 
nimibers,  the  time  to  process  the  grid  increases  for  the  FOUR-POINT  and 
EIGHT-POINT  methods. 
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Figure  2.1. Grids  showing  how  three  of  the  gridding  algorithms  interpolate 
a  cell  value. 
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Determining  Cell  Size  for  the  Grid 

Prior  to  the  gridding  process,  you  should  determine  an  appropriate  grid  cell 
size.  A  large  cell  size  (few  rows  and  columns  in  the  grid)  is  less  time  consum- 
ing but  may  be  also  less  accurate  than  a  small  cell  size.  You  can  use  a  large 
cell  size  for  a  first  run.  The  TESTGRID  command  allows  you  to  visually 
determine  what  cell  size  would  be  most  useful.  A  very  general  rule  is  that 
there  should  be  no  more  than  1  point  to  a  cell  and  points  not  more  than  10 
cells  away  from  each  other. 

Setting  the  "Roving  Window" 

Along  with  the  cell  size,  the  FOUR-POINT  and  EIGHT-POINT  gridding 
algorithms  will  ask  for  the  size  of  a  moving  window,  called  a  "roving  win- 
dow" in  the  MOSS  manual.  This  is  the  number  of  cells  on  all  sides  of  the 
target  cell  that  the  program  will  search  for  other  point  values  to  use  in  the 
interpolations.  A  very  general  rule  is  to  use  one-half  the  number  of  rows  in 
the  grid  pattem  for  the  size  of  the  roving  window. 


CHAPTER  TWO:  OSS  Applications  lor  Coal  Ro*oureo  Analyso*       41 


uajkayagaiaaaMMMi  jBMMMi.iiM  sji.miim 


TABLE  2.1.  Comparison  off  the  Gridding  Algorithms 


ALGORITHM 


ADVANTAGES 


DISADVANTAGES 


FOUR-POINT  Works  best  for  clustered  data 


Slow 


HGHT-POINT 


Only  method  that  allows  use  of  faults 
Produces  smooth  contour  maps 


Slow 


KRIGING 


Statistically  sound 


Difficult  to  learn 


QuiNrnc 

SPLINE 


Fastest 


Increasing  distortion  outside 
of  control  points.  Sees  trends 
where  data  is  sparse.  Often  has 
values  that  exceed  minimum/ 
maximum  values  of  the  data. 


POLYGON 


CTOG 


You  can  contour  by  hand  and 

turn  your  map  into  a  grid  for  analyses 


Can't  be  used  for  contour  maps 
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FLOW  CHART  2J3.  Generate  a  coal  thickness  grid 


•  VARIOGRAM  to  identify 
dupiicate  points 


T 


•  SELECT  points,  exciuding 

dupiicate  points 

•  SAVE  for  new  map 


J 


•  TESTGRiD  to  determine 
optimum  ceii  size 


I 


•  SELECT  by  thicioiess, 
exciude  zero  values 


I 


GRID  ttiicicness  vaiues 


PLOT  grid 
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COAL  EXERCISE  3.  Generate  a  coal  thickness  grid 

a.  Check  for  duplicate  points: 

(1)  FREE  all  maps  in  the  active  table  except  the  boundary  map 
(SJ.BORDER),  which  was  generated  in  Exercise  2.  WINDOW  to 
SJ.BORDER. 

(2)  SELECT  the  MYCOAL.PTS  map. 

(3)  Use  VARIOGRAM  on  the  Active  ED  for  MYCOAL.PTS  to  see  if 
there  are  any  duplicate  points  in  the  MYCOAL.PTS  map.  Use  all 
the  defaults.  Use  only  one  variogram  direction  with  angle  pairs  of 
22  and  0. 


In  this  case,  no  duplicate  points  are  found.  If  there  wsre  duplicate  points,  the  program  would 
have  listed  them  by  the  item  numbers.  You  would  SELECT  by  ITEM,  using  the  backslash  to 
"not  select"  one  of  each  pair  of  duplicate  points.  You  would  then  have  to  SAVE  this  map  and 
do  furtiier  analyses  with  the  new  map. 

h.  Use  TESTGRID  to  display  a  grid  over  the  study  area.  Use  meters  instead 
of  acres.  PLOT  the  MYCOAL.PTS  map  over  the  displayed  grid  to  check  iT 
the  cell  size  is  appropriate  for  the  gridding  process. 

c.  Generate  the  grid: 

(1)  BSEARCH  the  MYCOAL.PTS  map  lo  obtain  thickness  values  that  are 
greater  than  zero. 


Thickness  of  all  coal  beds  are  already  totaled  for  the  attribute  COAL.  THICK.  If  they  were  not 
tolled,  you  could  either  (a)  grid  each  bed  separately  and  then  use  the  MATH  command  to  add 
the  grids;  or  (b)  use  th&€OMPUTE  command  to  add  the  coal  thickness  values  for  the  beds 
and  create  a  new  attribute  for  the  total  sum. 

(2)  GRID  the  thickness  values  using  the  EIGHT-POINT  and  the  QUIN- 
TIC  SPLINE  algorithms. 

*  Use  a  250-meter  cell  size  'ir  all  map  analyses  in  these  exercises. 

*  For  the  EIGHT-POINT  gna,  use  a  roving  window  that  is  one-half  th. 

number  of  rows,-  for  the  Quintic  Spline,  use  a  window  of  7. 

(3)  DESCRIBE  the  grids  by  header  to  determine  if  they  have  cell  values 
that  are  not  within  the  minimum/maximum  values  of  the  data.  In 
MAPS,  RENUMBER  values  below  the  minimum  and  above  the 
maximimi  values. 
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(4)  PLOT  the  grids  and  graphically  compare  the  EIGHT-POINT  grid  with 
the  Quintic  Spline  grid. 
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FLOW  CHART  2A.  Evaluate  a  thickness  grid 


•  PLOT  data  points 

•  Label  thickness  values 


T 


•  PLOT  grid  and  look  for 
tight  concentric 
contours 


T 


•  PROFILE  and  look  at 
cross-sections 


I 


Display  grid  in  3D 
to  find  spikes 
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COAL  EXERCISE  4.  Evaluate  a  thicicness  grid. 

How  do  you  know  if  the  thickness  grid  you  calculated  is  correct? 

For  some  users,  this  "correctness"  is  a  limit  on  the  difference  between  the 
observed  and  calculated  thickness  values  in  the  grid.  For  other  users,  "cor- 
rectness" is  the  best  grid  one  can  obtain  given  certain  Umitations,  such  as 
CPU  time  or  the  amount  of  storage  available  (i.e.,  the  size  of  the  output  grid). 
The  geologist  must  also  determine  if  the  grid  is  geologically  acceptable. 

a.  Here  are  some  techniques  available  that  help  evaluate  the  output  grid: 

(1)  PLOT  the  data  points  and  label  (LEGEND)  them  by  the  thickness 
attribute.  See  if  any  thickness  values  look  out  of  line  with  neighbor- 
ing values. 

(2)  You  can  spot  data  errors  when  the  grid  is  plotted  and  there  are  tightly 
concentric  contours  around  a  point.  In  MOSS,  PLOT  the  EIGHT- 
POINT  grid  to  check  for  this. 

(3)  Use  the  PROFILE  command  to  look  at  the  grid  in  several  directions  in 
cross  section.  Do  the  cross  sections  look  acceptable? 

(4)  You  can  also  spot  data  errors  as  spikes  (positive  or  negative)  in  a  3-D 
display.  In  MAPS,  use  the  3D  command  to  display  the  EIGHT-POINT 
grid  to  locate  any  spikes. 
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FLOW  CHART  2.5.  Customize  a  thickness  grid 


Options 
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COAL  EXERCISE  5.  Customize  a  thickness  grid 

a.  Assuming  item  1  in  the  point  map  is  an  erroneous  ciata  point,  correct  or 

delete  it  by  one  of  the  following  methods  and  then  regrid  the  thickness: 

METHOD  A 

( 1 )  SELECT  the  point  map  (MYCOAL.PTS)  by  item,  using  the  backslash 
to  "not  select"  the  item  1. 

(2)  SAVE  this  map  and  give  it  a  new  name. 

(3)  SELECT  saved  map  and  reGRID  the  thickness,  using  the  EIGHT- 
POINT  method 

METHOD  B 

(1)  Use  EDIT  ATT  on  the  point  map  (MYCOAL.PTS)  to  edit  a  thickness 
value  for  item  1,  giving  it  a  new  value  of  1.0. 

(2)  ReGRID  the  edited  map,  using  the  EIGHT-POINT  method. 

b.  DESCRIBE  the  grid  by  header  to  see  if  there  are  negative  values.  If  so, 
RENUMBER  the  grid  so  that  all  the  negative  values  are  0  (zero). 

c.  After  you  have  completed  the  customization  of  the  thickness  grid,  rename 

it  SJTOK.GRD;  we  will  use  it  in  the  remaining  exercises  in  this  manual. 

d.  PLOT  the  grid  (SJTHK.GRD),  data  points  (MYCOAL.PTS),  and  lease  map 
(SJ.LEASE)  on  the  screen. 

e.  Clean  up  files: 

( 1 )  Keep  your  final  EIGHT-POINT  grid  (SJTHK.GRD). 

(2)  DELETE  other  thickness  grids  from  your  working  area. 

(3)  FREE  unnecessary  maps  from  your  active  table. 
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FLOW  CHART  2.6.  Calculate  Reserve  Base  tonnage 
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WINDOW  to  thickness  grid 
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Lease  boundary  metliod 
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•  Intersect  lease  and 
sections 
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•  Rasterize  the 
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TABLE  2.2.  Examples  of  reports  from  the  TOTAL  command  for  the 
lease  and  for  the  sections  within  the  lease. 


TOTAL  REPORT 

MAPNAME 
SITHK.GRD 
BYPC.LEASE 


PAGE  1 


FACTOR      CELL  SIZE 
LOOO  15.444 


SUBJECTS 

TOTAL 

FREQUENCY 

yiEA 

PRODUCT 

AVERAGE 

LEASE  1234 

14658. 

650. 

10039. 

226377.149 

22.551 

TOTAL  14658.  650. 

(BACKGROUND ACRES   -     34718. 


10039. 


226377.149 


22.551 


TOTAL  REPORT 

MAPNAME 
SrrHK.GRD 
BYPC.SECnON 

SUBJECTS 


TOTAL 


FACTOR 
1.000 


CELL  SIZE 

15.444 


TOTAL 


FREQUENCY      AREA 


14658. 


:  BACKGROUND  ACRES 


34718. 


650. 


10039. 


PAGE  1 


PRODUCT     AVERAGE 


01;T23N.R13W 

1072. 

49. 

757. 

16558.980 

21.881 

02;T23N.R13W 

1044. 

42. 

649. 

16121.633 

24.854 

04;T23N.R12W 

1313. 

42. 

649. 

20274.015 

31.256 

05;T23N.R12W 

1749. 

40. 

618. 

27006.316 

43.716 

06;T23N.R12W 

1857. 

41. 

633. 

28678.471 

45.291 

07;T23N.R12W 

1089. 

36. 

556. 

16816.888 

30.247 

08;T23N.R12W 

1050. 

43. 

664. 

16220.969 

24.426 

09)T23N.R12W 

1008. 

46. 

710. 

15562.364 

21.906 

11;T23N.R13W 

618. 

36. 

556. 

9538.365 

17.156 

12;T23N.R13W 

606. 

42. 

649. 

9358.956 

14.428 

13,T23N.R13W 

746. 

49. 

757. 

11516.142 

15.218 

14,T23N.R13W 

599. 

42. 

649. 

9247.675 

14.257 

15;T23N.R13W 

783. 

49. 

757. 

12096.794 

15.985 

17;T23N.R12W 

469. 

48. 

741. 

7236.752 

9.762 

18;T23N.R12W 

610. 

42. 

649. 

9424.286 

14.529 

226377.149 


22.551 
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FLOW  CHART  2.7.  Generate  overburden,  stripping  mining-ratio 
and  coal-potential  grids 
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mir 


Generate  mining-ratio  grid 


Generate  overburden  grid 


J 


Multiple  Recovery  Factors?- 


1 

NO 

J 


Generate  mining-ratio  grid 


MATH:  (Overburden  grid  *  (1613.33  / 
Coal  Factor) )  divided  by  (thickness 
grid  times  factor  grid) 


•  MATH:  (Overburden  grid  *  (1613.33  / 
Coal  Factor) )  divided  by  (thickness 
grid  times  Recovery  Factor) 


T 


Generate  coal-potential  grid 


Group  mining-ratio 
values  for  high, 
medium,  low  potential 
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COAL  EXERCISE  6.  Calculate  the  Reserve  Base  tonnage 

After  you  have  edited  the  thickness  grid  and  the  grid  meets  your  require- 
ments, always  window  to  this  grid  for  the  remaining  processes.  Do  this 
because  it: 

•  provides  consistency  for  the  project 

•  is  time-saving  to  always  window  to  one  map 

•  meets  window  requirements  for  raster  analyses,  i.e.,  all  analyses  must 

be  with  maps  containing  the  same  number  of  rows  and  columns 

In  this  study  area,  the  coal  thickness  is  greater  than  5  feet  within  the  lease 
boimdary.  If  there  were  coal  less  than  5  feet  and  you  wanted  to  limit  the 
toimage  calculations  to  a  coal  thickness  greater  than  this  value,  you  would 
first  RENUMBER  the  thickness  grid,  assigning  the  value  of  0  to  all  cells 
values  below  4.99. 

a.  Calculate  the  Reserve  Base  tonnage  within  the  lease  boimdary: 

(1)  In  MOSS,  WINDOW  to  the  coal  thickness  grid  (SJTHK.GRD). 

(2)  SELECT  and  POLYCELL  the  lease  map  (SJ.LEASE)  using  the  same  cell 
size  as  the  thickness  grid  (SJTHK.GRD). 

•  Name  the  new  map  PC.LEASE 

(3)  DESCRIBE  the  thickness  grid  (SJTHK.GRD)  and  the  polycelled 
leasemap  (PC.LEASE)  to  make  sure  the  number  of  rows  and  columns 
are  the  same. 

(4)  In  MAPS,  WINDOW  to  the  thickness  grid  (SJTHK.GRD). 

(5)  Use  the  PAGE  command  to  increase  page  length  to  a  large  number 
(e.g.,  3000)  to  prevent  multiple  headings  in  the  tonnage  report  gener- 
ated in  the  next  step. 

(6)  Use  the  TOTAL  command  to  determine  the  average  thickness  for  the 
lease  area: 

•  Use  the  thickness  grid  (SJTHK.GRD) 

•  Use  a  factor  of  1 

See  Table  2.2  for  an  example  of  a  report  generated  by  the  TOTAL  command. 
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(7)  Calculate  the  Reserve  Base  (RB)  tonnage  for  the  lease  by  manually 
(with  hand  calculator)  multiplying  the  following  values: 

•  the  average  thickness  from  the  tonnage  report 

•  the  Coal  Factor  (1770  tons/acre-feet  of  coal) 

•  the  area  in  acres  obtained  by  using  the  AREA  command  in  MOSS 
on  the  lease  map  (SJ.LEASE) 

iir 

The  reason  we  use  the  acreage  from  the  vector  (MOSS)  side  ofGIS  is  thatrasterizing  (changing 
a  map  from  vector  to  cells)  often  changes  the  total  acreage  within  a  polygon  depending  on  the 
cell  si2B  used  A  cell  must  contain  over  50  percent  of  the  polygon  to  tje  included  in  the  cell 
map.  If  a  cell  contains  less  than  50  percent,  it  is  considered  background 

b.  Calculate  the  Reserve  Base  tonnage  within  each  section  within  the  lease 
botmdary: 

(1)  In  MOSS,  SELECT  the  township/range/section  map  (LANDLINES). 
OVERLAY  the  lease  map  (SJ.LEASE)  and  the  landlines  map. 

•  Use  zero  characters  for  the  lease  and  16  characters  for  the 

landlines. 

•  Use  the  intersection  option  (menu  number  1). 

•  Name  the  output  map  OV.SECTION  (OV  for  overlay). 

(2)  SELECT  the  OV.SECTION  map.  POLYCELL  ..le  intersected  map, 
using  the  same  window  and  cell  size  as  the  coal  thickness  grid. 

•  Name  the  output  map  PC.SECTION. 

(3)  DESCRIBE  the  thicLness  grid  and  the  rasterized  section  map  to  make 
sure  the  number  of  rows  and  columns  are  the  same. 

(4)  Use  the  TOTAL  command  to  determine  the  average  thickness  for  the 
lease  area: 

•  Use  the  thickness  grid 

•  Use  a  factor  of  1 

•  Name  the  output  report  SECTION.RPT  (see  Table  2.3) 

c.  OPTIONAL:  Manually  calculate  the  Reserve  Base  toimage  for  each  section 
within  the  lease  bovmdary.  Use  a  Coal  Factor  of  1770,  the  average  th;  k- 
ness  as  shown  in  SECnON.RPT,  and  vector  acreage  for  the  sections. 
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COAL  EXERCISE  7.  Generate  overburden,  stripping  mining-ratio, 

and  coal  potential  grids 

In  the  SJCOAL  study  area,  a  DEM  is  not  available  to  use  in  generating  an 
overburden  grid.  However,  the  multiple  attributes  attached  to  the  data 
points  provide  overburden  and  interburden  values. 

a.  Generate  an  overburden  grid: 

(1)  First  create  a  new  attribute  for  the  sum  of  interbiurden  and  overburden: 

(a)  Using  the  COMPUTE  command  on  the  Active  ID  of  your  coal 
data  map  (MYCOALPTS),  combine  the  thickness  of  the 
overburden  and  the  thickness  of  the  interburden  to  create  a  new 
attribute  named  OVERINTB. 

(2)  Generate  the  overburden  grid  using  the  new  attribute: 

(a)  SELECT  and  GRID  the  new  attribute  OVERINTB. 

•  Window  to  the  thickness  grid  (SJTHK.GRD) 

•  Name  the  output  grid  OVERB.GRD 

•  Use  the  EIGHT-POINT  gridding  algorithm 

(3)  Display  the  overburden  grid  on  the  screen: 

(a)  FLOOD  areas  where  the  coal  bed  might  be  outcropping  (where 
overburden  is  zero).  There  may  not  be  outcrops  in  this  study  area. 

(b)  Areas  where  the  coal  beds  are  overlain  by  200  feet  or  less  of 
overburden  are  considered  to  have  potential  for  strip  mining.  To 
visually  determine  if  the  lease  has  this  potential,  FLOOD  the 
overburden  grid  (OVERB.GRD)  for  three  categories: 

5  through  60  feet 

60.00001  through  120  feet 

120.00001  through  200  feet 

PLOT  and  SE/  DE  the  rasterized  lease  over  the  overburder  display. 

b.  In  MAPS,  generate  a  stripping  mining-ratio  grid  using  a  Recovery  Factor 
and  a  Coal  Factor: 

(1)  As  mining  ratios  are  based  on  cubic  yards  of  overburden  per  ton  of 
:  recoverable  coal,  the  following  conversion  formula  was  calculated: 
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Cubic  Yards  of  rock  /  ton  of  coal  = 

43,560  cu.ft.  /  acre-feet  = 

27  cubic  ft.  /  cubic  yard 

1 6 13.33  cubic  yards/acre-feet 

(2)  Mining  ratios  are  calculated  using  the  following  formula: 

MR  =  To  M  161333 /CF) 

( Tc  )  *  RF 

where  MR  =  Mining  Ratio 

To     =  thickness  of  overburden 

1613.33     =  cubic  yards/acre-feet 

CF     =  Coal  Factor 

Tc     =  thickness  of  coal 

RF     =  Recovery  Factor 

Use  the  IVIATH  command  and  substitute  grids  in  the  formula.  Use  a 
Coal  Factor  of  1 770  (for  subbituminous  coal)  and  a  Recovery  Factor  of 
.85.  Name  the  new  grid  RATIO.GRD: 

MATH  (  OVERB.GRD  *  (  1613.33  /  1770  )  )  /  , 
(  SJTHK.GRD   *  .85  )  FOR  RATIO.GRD 

c.  Look  at  the  mining-ratio  grid  graphically: 

(1)  In  MAPS: 

(a)  FLOOD  the  ratio  grid  using  different  colors  to  display  the  follow- 
ing four  categories: 

0  through  5 

5.00001  through  10 

10.00001  through  15 

Greater  than  15.00001 


For  these  exercises  we  have  used  3  large  cell  si2s  in  order  to  conserve  computing  time. 
With  a  smaller  cell  size,  the  step-like  interfyces  would  tx  less  noticeable 
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(b)  PLOT  the  rasterized  lease  map  (PC.LEASE)  over  the 
mining-ratio  display. 

(c)  In  MOSS,  SELECT  and  PLOT  the  mining-ratio  grid 
(RATIO.GRD),  using  a  contour  interval  of  five. 

d.  Generate  a  mining-ratio  grid  using  multiple  Recovery  Factors: 

•  Recovery  Factor  is  .85  v/here  overburden  is  equal  or  less  than  100  feet 

•  Recovery  Factor  is  .50  where  overburden  is  greater  than  100  feet 

( 1 )  RENUMBER  the  overburden  grid  (OVERB.GRD)  to  reflect  the  above 
parameters.  Name  the  new  grid  FACTOR. GRD. 

OPTIONAL:  Validate  the  output  by  flooding  the  overburden  grid 
(OVERB.GRD),  assigning  a  color  from  0  feet  through  100  feet.  Then 
shade  the  factor  grid  (FACTOR.GRD),  assigning  1  to  .85.  The  area 
covered  by  these  two  should  be  the  same. 

(2)  Use  the  MATH  command  and  substitute  grids  in  the  mining-ratio 
formula.  Use  a  Coal  Factor  of  1770  (for  subbituminous  coal)  and  the 
Recovery  Factor  grid  (FACTOR.GRD).  Name  the  new  grid 
RATIOGRD.2: 

MATH    (    OVERS. G!RD   *    (    1613.33    /    1770    ))     /    , 
(    SJTHK.GRD    *    FACTOR.GRD    )    FOR   RATIOGRD.2 

OPTIONAL:  VaUdate  the  output  by  first  flooding  the  factor  grid 
(FACTOR.GRD),  assigning  a  color  to  .85.  All  cells  within  this  area  should 
have  the  same  mining-ratio  value  for  the  RATIO.GRD  and  RATIOGRD.2. 
Query  several  cells  at  random,  especially  at  the  interface  between  .85  and  .5 
and  see  if  the  cells  for  both  grids  have  the  same  mining-ratio  value. 

(3)  FLOOD  the  new  mining-ratio  grid  (RATIOGRD.2),  using  different 
colors  to  display  the  same  f,  ;x  categories  as  you  did  for  the  first  ratio 
grid  (RATIO.GRD).  PLOT  and  SHADE  the  rasterized  lease  grid 
(PC.LEASE)  over  this  ratio  grid. 

e.  Generate  a  coal  development  potential  grid: 

The  definition  of  coal  development  potential  changes;  it  is  based  on 
present-day  economic  and  technological  criteria  as  well  as  geologic  crii.  j 
withinthe  study  area.  Use  the  following  definitions  for  this  exercise: 
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MINING  RATIO DEVELOPMENT  POTENTIAL 

0  to  5:1  High 

5:1  to  15:1  Medium 

Greater  than  15:1  Low 

(e.g.,  5: 1  meaning  5  cubic  yards  of  rock  per  1  ton  of  coal) 

(1)  RENU2VIBER  the  cell  values  in  the  mining-ratio  grid  (RATIO.GRD) 
and  group  them  into  the  three  potential  categories: 

•  Assign  cell  value  1  to  high  potential,  2  to  medium,  and  3  to  low. 

•  Name  the  output  grid  COALPOT.GRD. 

(2)  FLOOD  the  coal  potential  grid  (COALPOT.GRD)  with  different  colors. 
The  cell  values  are  1, 2,  and  3.  SHADE  the  rasterized  lease  map 
(PC.LEASE)  over  the  grid. 
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FLOW  CHART  Z^S.  Generate  Surface-Minable  Coal  (SMC)  grid  and 

calculate  SMC  tonnage 


Stripping  mining-ratio  grid 
(single  or  multiple 
Recovery  Factors) 


Generate  GoaNlmitatlon  grids 


'  Grid  analytical  values 
(asii,  sulfur,  etc.) 
—  —  -and  —  —  — 

>  Digitize  and  rasterize  oUier 
limitations  sucii  as  bum  or 
or  mined  areas  or  cliannels 


T 


Generate  minable-coal  grid 


'  Cross  grids,  assign  values 
1  to  cells  witliln  limits 


f 


Generate  Surface  IMinabie-Coal  (SMC)  grid 


•  Cross  ratio  grid  with  minable-coal 
grid,  assign  value  1  to  cells 
within  limits 


T 


Generate  SMC  grid  contained  within  lease 


'  Cross  SMC  grid  with  rasterized  lease 


I 


•  Label  SMC-lease  grid:  maintain  cell 
values;  enter  subject  names 


I 


•  Total  thickness  grid  by  SMC-lease 
grid  for  tonnage  report 


aot  covered  la  this  training  manual 
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COAL  EXERCISE  8.  Generate  a  Surface-Minabie  Coal  (SMC)  grid  and 

calculate  SMC  tonnage 

A  grid  of  suiface-minable  coal  (SMC)  can  be  generated  to  contain  only  those 
grid  cells  which: 

(a)  meet  specific  analytical  criteria  and  other  coal  limitations  (mined 
or  burned  areas,  chaimels); 

(b)  meet  a  specific  mining-ratio  limit;  and 

(c)  are  within  the  lease  boundary. 

The  resultant  grid  can  then  be  used  with  the  coal  thickness  grid  to  calculate 
the  coal  tonnage  of  surface-minable  coal  within  the  lease. 

a.  Generate  a  minable-coal  grid  that  meets  specific  analytical  criteria  and 
other  coal  limitations: 

(1)  Generate  grids  for  analytical  factors  that  limit  the  coal  that  can  be 
mined.  Use  the  SULFUR  and  ASH  attributes  in  the  map 
COAL.SAIVIP: 

•  WINDOW  to  the  thickness  grid  (SJTHK.GRD). 

•  Use  the  Quintic  Spline  option  since  there  are  only  nine  points  in 

the  map. 

•  Use  the  same  cell  size  as  the  thickness  grid 

•  RENUJVIBER  grid  if  there  are  negative  values. 

•  Name  the  sulfur  grid  SULFUR.GRD. 

•  Name  the  ash  grid  ASH.GRD. 

(2)  CROSS  the  ash  and  sulfur  grids  so  that  cell  values  meeting  thecriteria 
have  a  value  of  1,  and  ceil  values  above  the  criteria  have  a  value  of 
zero. 

•  Use  values  of  0  through  0.6  percent  as  an  acceptable  range  for 

sulfur  content. 

•  Use  values  of  0  through  20  percent  as  an  acceptable  range  for 

ash  content. 

•  Name  the  output  grid  MINABLE  for  minable  coal. 
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EXAMPLE: 

CROSS  sulfur.grd  WITH  ash.grd  FOR  minable , 

ASSIGN  1  TO  0  THRU  .6  AND  0  THRU  20, 

ASSIGN  0  TO  .60001  THRU  9999  AND  20.00001  TH  9999 


The  first  assigned  values  are  for  the  sulfur  grid  (SULFUR.  GRD). 
The  second  assigned  values  afB  for  the  ash  grid  (ASH.GRD).  The  "AND"  between  these  two 
phrases  is  a  Boolean  option  and  requires  that  the  first  condition  must  tie  met  AND  the  second 
condition  must  be  met  for  the  new  grid,  IvIINABLE. 

(3)  FLOOD  the  minable-coal  grid  (MINABLE),  assigning  a  color  num- 
ber to  cells  with  a  value  of  1. 

b.  Generate  Surface-Minable  Coal  (SMC)  grid: 

(1)  CROSS  the  minable-coal  grid  (MINABLE)  with  the  mining  ratio  grid 
(RATIO.GRD): 

•  Assign  value  of  1  for. 

1  through  1  cell  values  as  the  acceptable  range  for  the  minable 
coal. 

0.5  through  5  cell  values  as  an  acceptable  range  for  the  mining 
ratio. 

•  Name  the  output  grid  SMCOAL  for  surface-minable  coal. 

(2)  FLOOD  the  minable  coal  (MINABLE),  assigning  a  color  to  cells  with  a 
value  of  1.  SHADE  the  surface-minable  coal  (SMCOAL),  assigning  a 
shade  pattern~to  cells  with  a  value  of  1. 

c.  Find  the  Surface-Minable  Coal  within  the  lease  boundary: 

(1)  CROSS  the  rasterized  lease  (PC.LEASE)  with  the  Surface-Minable  Coal 
grid  (SMCOAL): 

•  Assign  a  value  of  1  for: 

1  through  1  cell  values  in  the  lease  map. 

1  through  1  cell  values  in  the  Surface-Minable  Coal  grid. 

•  Name  the  new  grid  SMCOAL.LSE  (LSE  meaning  lease). 
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(2)  This  new  grid  is  without  a  label  (subject)  and  will  result  in  an  error 
message  when  using  the  TOTAL  command.  Use  the  LABEL  com- 
mand to  add  a  subject: 

•  Use  16  characters  for  the  subject  (e.g.,  'SMC  in  lease'). 

•  Value  will  remain  "I." 

(3)  Look  at  the  grids  on  the  screen: 

FLOOD  the  minable  coal  (MINABLE),  assigning  a  color  to  cells  with 
a  value  of  1. 

SHADE  the  surface-minable  coal  (SMCOAL),  assigning  a  shade 
pattern  to  cells  with  a  value  of  1. 

SHADE  the  surface-minable  coal  within  the  lease  boundary 
(SMCOAL.LSE),  assigning  another  shade  pattern  to  cells  with  a 
value  of  1. 

PLOT  the  rasterized  lease  boimdary  (PC.LEASE)  in  a  contrasting 
color. 

CONTOUR  the  thickness  grid  (SJTHK.GRD)  with  an  interval  of  2  — 
and  with  a  contrasting  color. 

d.  Calculate  the  toimage  of  surface-minable  coal  using  the  TOTAL  com- 
mand: 

•  Use  the  coal  thickness  grid  (  SJTHK.GRD)  and  the  surface  minable 

grid  within  the  lease  (SMCOAL.LSE). 

•  Use  a  factor  of  1  to  get  the  average  thickness  (see  Table  2.3). 

•  Repeat  the  command  using  a  coal  factor  of  1 770  tons/acre-  feet  of 

coal  (see  Table  2.3). 
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TABLE  2^.  Example  off  report  from  the  TOTAL  command  ffor  the 
surface-minabie  coal  within  tiie  lease  boundary 


TOTAL  SITHK-GRD  BY  SMCOAL.LSE  LABEL  TONS  FACTOR  1 


TOTAL  REPORT 

MAPNAME 
SJTHKGRD 
BY  SMCOAL.LSE 


PAGE  1 


FACTOR      CELL  SIZE 
LOGO  15.444 


SUBJECTS 

TOTAL 

FREQUENCY 

AREA 

TONS 

AVERAGE 

SURF.MINCOALIN 

7452. 

261. 

4031. 

115086.003 

28.551 

TOTAL  7452  261. 

(BACKGROUND ACRES   -     40726. 1 


4031. 


115086.003       28.551 


TOTAL  SJTHK.GRD  BY  SMCOAL.LSE  LABEL  TONS  FACTOR  1 770 


TOTAL  REPORT 

MAPNAME 

SJTHKGRD 

BY  SMCOAL.LSE 


FACTOR       CELL  SIZE 
1770.000        15.444 


PAGE  1 


SUBJECTS 

TOTAL 

FREQUENCY 

AREA 

TONS 

AVERAGE 

SURF.MINCOALIN 

7452. 

261. 

4031. 

203702226.089 

50535.091 

TOTAL  7452.  261. 

( BACKGROUND  ACRES   =     40726.) 


4031. 


203702226.089        50535.091 
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FLOW  CHART  2.9.  Prepare  the  DEM  for  map  analyses 


MAGNET  to  upload 
magnetic  tape  Into 
your  directory 


J 


IMPORT  DEM  Into  MAPS 


J 


RePROJECT  DEM  if  in 
latitlude/iongitude 


I 


MERGE  OEMs  If  required 


I 


•  CUT  DEM  by  tliici(ness  grid  window 

•  RESAMPLE  DEM  to  ceil  size  Of  thickness  grid 

•  MATH  to  convert  DEM  from  meters  to  feet,  using 

INTeger  function 

•  LABEL  converted  DEM  witli  correct  description 


not  covered  In  this  training  manual 
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Loading  the  DEMs  for  Use  in  MOSS/MAPS 

Loading  the  DEM  tapes  and  importing  the  files  into  MAPS  is  not  covered  in 
this  manual.  However,  here  is  a  brief  description  of  the  steps  you  must  take. 

Upload  the  DEMs: 

Inform  the  computer  operator  in  your  office  which  DEM  quadrangle  you 
want  uploaded,  e.g.,  B3444106.  Usually  the  computer  operator  will  upload 
the  DEM  file  from  a  tape  into  your  directory.  If  this  is  not  the  case,  then  take 
the  following  steps: 

Login  into  PRIME  and  assign  the  correct  tape  drive  number.  Then  use  the 
MAGNET  program  to  transfer  the  data  from  the  tape  into  a  file  in  your 
directory. 

Example: 

ASSIGN  MTO  [or  MTl,  depending  on  which  tape  drive  is  being  used] 

MAGNET 

Read 

MTU  No.:  0   [should  be  same  drive  nvunber  used  above] 

File  No.:  1   [or  the  position  of  the  file  on  the  tape] 

Logical  Record  Length:  1024 

Blocking  Factor:  4 

Ascii,  binary,  etc:  ASCII 

Diskfile:  dem.name  [enter  the  output  filename] 

QUIT     [wait  imtil  a  '>'  prompt  appears  before  you  exit  MAGNET] 
UNASSIGN  MTO  [rewinds  the  tape] 

Importing  the  DEM  Data  into  MAPS 

Use  the  IMPORT  command  in  MAPS  to  generate  a  map  from  the  C  ir  .  file. 
The  FORMA        ould  be  "DEM  7.5 "  and  the  map  TYPE  should  be    8. " 

RePRO  JECT  ana  2V1ERGE  several  DEM  maps  together,  if  required,  for  your 
study  area. 
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Given  Data  for  the  DHCOAL  Study  Area 

In  the  DHCOAL  study  area,  there  are  45  drill  holes.  The  DEM  is  in  meters, 
has  a  30-meter  cell  size,  and  covers  a  larger  area  than  the  drill  holes.  Thus, 
the  DEM  must  be  converted  to  feet,  resampled  to  a  larger  cell  size,  and  cut  to 
fit  the  size  of  the  study  area. 

Master  maps  in  the  MOSS  project  DHCOAL  are: 


Master  Map  Name        Map  Type        Map  Theme 


DHPLS 

3  (polygon 

DH.PTS 

1  (point) 

DHDEMM 

8  (ceU) 

DHSTR.GRD 

8  (ceU) 

DHTHK.GRD 

8(ceU) 

Township,  range,  section 

Coal  data  for  45  drill  holes 

DEM  in  meters 

Structure  grid  for  the  top  of  a  coal  bed 

Thickness  grid  for  the  coal  bed 
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EXERCISE  9.  Prepare  the  DEM  for  map  analyses 

For  the  following  exercises,  you  need  to  change  projects  within  the  Master 
directory  ZVIINTRAIN.  OPEN  to  the  project  DHCOAL. 

a.  In  MOSS,  DESCRIBE  the  DEM  (DHDEM.M)  by  header.  Note  that 
the  minimum  and  maximum  values  are  1864  meters  and  3064  meters. 
DESCRIBE  the  DEM  by  projection.  Note  that  the  map  is  in  the  UTM 
projection. 

b.  WINDOW  to  the  thickness  grid  (DHTHK.GRD).  CUT  the  DEM  by  the 
window  so  that  it  has  the  same  number  of  rows  and  columns  as  the 
thickness  grid: 

•  Name  the  cut  map  DEM.CUT. 

c.  RESA2VIFLE  the  cut  DEM  so  that  it  has  the  same  cell  size  as  the  thickness 
grid: 

•  Name  the  resampled  map  DEM.RES. 

d.  Convert  the  resultant  grid  from  meter  to  feet  using  the  MATH  command: 

•  Use  the  INT  function. 

•  Multiply  by  a  conversion  factor  of  3.28083989. 

•  Name  the  converted  map  DEM.FT. 

e.  DESCRIBE  (by  header)  the  converted  DEM  (DEM.IT)  and  the  thickness 
grid  (DHTHK.GRD)  to  make  sure  the  number  of  rows  and  columns  are  the 
same. 

f .  Use  the  LABEL  command  in  MAPS  to  change  the  description  of  the 
converted  DEM  so  that  it  reads  "DEM  in  feet." 

g.  Clean  up  the  files: 

(1)  Keep  the  DEM,  in  feet  (DEM.FT). 

(2)  DELETE  the  following  grids: 

DEM.CUT 
DEM.RES 
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FLOW  CHART  2.10.  Compare  collar  elevations  to  DEM  elevations 


•  WINDOW  to  DEM,  In  feet 


I 


Compare  collar  elevation  with  DEM  elevation: 

•  Rasterize  collar  elevation 

•  Rasterize  subjects 

•  Subtract  collar  elev.  from  DEM 

•  Categorize  resultant  cell  values 

•  increase  page  length 

•  Total  rasterlzed  subjects  with  categorized  map  to 

generate  report  showing  differences 
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COAL  EXERCISE  1 0.  Compare  collar  elevations  to  DEM  elevations 

You  will  now  use  the  following  sequence  of  commands  to  compare  the  collar 
elevations  of  the  drill  holes  to  the  DEM  elevations  at  the  same  point  location: 

POLYCELL  (rasterize)  observed  values 

POLYCELL  (rasterize)  subjects 

MATH  (subtract) 

CATEGORIZE 

PAGE 

TOTAL 

a.  Generate  the  rasterized  maps: 

(1)  In  MOSS,  WINDOW  to  the  DEM  in  feet  (DEM.Fr). 

(2)  SELECT  and  POLYCELL  the  point  map  (DH.PTS)  twice,  once  for  the 
attribute  option  and  once  for  the  subject  option. 

(a)  The  attribute  option: 

•  Name  the  output  grid  PC.COLLAR. 

•  Use  same  cell  size  as  thickness  grid  (DHTHK.GRD). 

•  Choose  option  2  (discrete  map  Type  7)  for  "new  cell  map 
type." 

NOTE:  This  must  be  a  Type  7. 

•  Choose  option  6  (multiple  attribute  option)  for  "cell  assign 
ment  option." 

•  Choose  the  attribute  COLLAR.ELV. 

(b)  The  subject  option: 

•  Name  the  output  grid  PC.SUBJ. 

•  Use  same  cell  size  as  thickness  grid  (DHTHK.GR 

•  Choose  option  2  (discrete  map  Type  7)  for  "new  ceil  m.  o 
type." 

-  —    -   -   Choose  option  2  (subject  option)  for  "cell  assignment  option." 
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b.  In  IVIAPS,  use  the  MATH  command  to  subtract  the  DEM  grid  (DEM.FT) 
from  the  collar  elevation  grid  (PC.COLLAR): 

•  Name  the  residual  grid  DIFF. 

c.  Obtain  a  listing  of  differences  in  the  elevations: 

( 1 )  CATEGORIZE  the  residual  grid  (DIFF)  to  create  a  Type  7  grid;  a  Type 
7  grid  is  needed  for  the  TOTAL  command: 

*  Name  the  output  grid  DIFF.T7  (T7  for  Type  7). 

(2)  Set  the  PAGE  length  to  a  large  number  (e.g.,  3000)  so  that  the  gener- 
ated report  will  be  continuous,  i.e.,  without  multiple  headings. 

(3)  TOTAL  the  categorized  map  (DIFF.T7)  with  the  polycelled  subject  grid 
(PC.SUBJ): 

•  Use  the  'LABEL'  option  for  the  report. 

(4)  In  PRIMOS,  SLIST  the  report,  which  lists  the  elevation  differences 
(column  marked  TOTAL')  and  the  subjects.  See  Table  2.5  for  an 
example  table.  The  differences  may  be  found  under  the  header  'TO- 
TAL.' 

OPTIONAL:  Find  only  those  differences  that  are  greater  than  plus/minus  50 
feet  in  elevation. 

(1)  CATEGORI2aE  the  residual  grid  (DIFF)  twice,  once  to  create  grid  for 
plus  values  and  again  for  minus  values: 

•  Use  the  FROM  option  (FROM  50  TH  9999). 

Name  the  output  grid  DIFF.T7P  (T7  for  Type  7  and  P  for 
'plus  values'). 

•  Use  the  FROM  option  (FROM  -9999  TH  -50) 

Name  the  output  grid  DIFF.T7M  (T7  for  Type  7  and  M  for 
'minus  values'). 

(2)  Set  the  PAGE  length  to  a  large  number  (e.g.,  3000). 

(3)  TOTAL  each  of  the  categorized  maps  (DIFF.T7P  and  DIFF.T7M)  with 
the  polycelled  subject  grid  (PC.SUBJ): 

•  Use  the  LABEL  option  to  generate  the  reports. 

d.  Clean  up  the  files: 
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( 1 )  DELETE  the  following  grids: 

PC.COLLAR  DIFF.T7 

PC.SUBJ  DIFF.T7M 

DIFF  DIFF.TTP 

(2)  FREE  unnecessary  maps/grids  from  your  active  table. 
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TABLE  2.4  Summary  table  for  the  differences,  in  feet,  of  collar  elevation 
and  the  DEM  elevation. 

NOTE:  The  difference  in  elevation  is  listed  under  the  header  'TOTAL', 


TOTAL  REPORT 

PAGE  1 

MAPNAME 

FACTOR 

CELL  SIZE 

Dll-i--2 

1.000 

15.444 

BYPC.SLFBJ 

SUBJECTS 

TOTAL 

FREQLnENCY 

AREA 

PRODUCT 

AVERAGE 

DH-1 

39. 

15. 

602.319 

39.000 

DH-10 

69. 

X, 

15. 

1065.642 

69.000 

DH-12 

33. 

i. 

15. 

509.655 

33.000 

DH-13 

17. 

IS. 

262.549 

17.000 

DH-14 

^. 

1* 

15. 

-61.776 

-4.000 

DH-15 

30. 

!■ 

15. 

463.323 

30.000 

DH-16 

-19. 

1, 

15. 

-293.438 

-19.000 

DH-18 

109. 

1* 

15. 

1683.405 

109.000 

DH-19 

33. 

1, 

15. 

509.655 

33.000 

DH-2 

89. 

Am 

15. 

1374.524 

89.000 

DH-20 

-134. 

1* 

15. 

-2069.508 

-134.000 

DH-21 

338. 

li 

15. 

5220.101 

338.000 

DH-22 

6. 

1, 

15. 

92.665 

6.000 

DH-23 

3. 

15. 

46.332 

3.000 

DH-24 

55. 

!■ 

15. 

849.425 

55.000 

DH-25 

-17, 

Xt 

15. 

-262.549 

-17.000 

DH-26 

177. 

1, 

15. 

2733.603 

177.000 

DH-27 

-183. 

J,, 

15. 

-2826.268 

-183.000 

DH-29 

■39. 

1, 

15. 

-602.319 

-39.000 

DH-3 

29. 

1, 

15. 

447.878 

-29.000 

DH-30 

-68. 

1, 

15. 

-1050.198 

-68.000 

DH-31 

108. 

1. 

15. 

1667.961 

108.000 

DH-32 

243. 

15 

3752.913 

243.000 

DH-33 

84. 

15. 

1297.303 

84.000 

DH-34 

193. 

15. 

2980.709 

193.000 

DH-35 

■36. 

15. 

-555.987 

-36.000 

DH-36 

-86. 

15. 

-1328.191 

-86.000 

DH-37 

33. 

15. 

509.655 

33.000 

DH-38 

-4. 

15. 

-61.776 

-4.000 

DH-39 

0. 

15. 

0.000 

0.000 

DH-4 

-158. 

15. 

-2440.166 

-158.000 

DH-40 

■27. 

15. 

-416.990 

-27.000 

DH-41 

■23. 

15. 

-355.214 

-23.000 

DH-42 

■29. 

15. 

-447.878 

-29.000 

DH-43 

3. 

15. 

46.332 

3.000 

DH-44 

36. 

15. 

555.987 

36.000 

DH-45 

29. 

15. 

447.878 

29.000 

DH-5 

158. 

15. 

2440.166 

158.000 

DH-6 

56. 

15. 

864.869 

56.000 

DH-7 

11. 

15. 

169.885 

11.000 

DH-8 

-3. 

15. 

-46.332 

-3.000 

DH-9 

-13. 

15. 

-200.773 

-13.000 
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FLOW  CHART  2.1 1.  Use  DEM  to  generate  an  overburden  grid. 


WINDOW  to  thickness  grid 


i 


•  Subtract  the  structure  grid  from  the  DEM 
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COAL  EXERCISE  1 1 :  Use  DEM  to  generate  an  overburden  grid 

a.  Generate  an  overburden  grid  using  the  structure  grid  (DHSTR.GRD) 
and  the  DEM  in  feet  (DEM  JT): 

(1)  In  MAPS,  WINDOW  to  the  thickness  grid  (DHTHK.GRD).* 

(2)  Use  the  MATH  command  to  subtract  the  structiire  grid 
(DHSTR.GRD)  from  the  DEM  (DEM-FT): 

•  Name  the  resultant  grid  OVERB.DEM. 

(3)  DESCRIBE  the  overburden  grid  (OVERB.DEM)  by  header  to  determine 
the  mLQimtrai  and  maximum  values.  Note  that  there  are  negative 
values,  indicating  that  there  is  an  outcrop  area.  RENUMBER  the 
negative  values  to  zero. 

(4)  Use  the  FLOOD  command  to  display  where  the  outcrop  area  is 
located. 

(5)  CONTOUR  the  renumbered  overburden  grid  over  the  outcrop  area 
you  flooded  in  the  previous  step. 
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Appendix 


Answers:  Coal  Exercise  1 

Step  A 


MOSS  MINTRAIN  usemane 
OPEN  SJCOAL 
FREE  ALL 

Stepb 

list  master 

Stepc 

UTILITY 

7 

SJCOAL 

1 

SYSTEM  SLIST  SJCOAL. BL 

Stepd 

describe  coal.pts  at 

4 
1 
SELECT  COAL.PTS  ALL 

WINDOW  #  /#  =  Active  table  ID  number] 

ERASE    -   PLOT  # 

QUERY 

/*   CROSSHAIRS 

0 

REPORT  # 

COAL.RPT 

Y 

N 

3 

BEDS 

3 

COAL. THICK 

3 

TOT.OVERB 

1 

SYSTEM  SLIST  COAL.RPT 
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Stepe 

bsearch  coal.pts 

tot.overb  lt  30  and  coal. thick  gt  15 

erase  -  plot  # 

legend  #  label 

4 

Stepf 

select  s j. lease  all 

SELECT  WSA  SUBJECT 

\NDLL 

SELECT  LANDLINES  ALL 

SELECT  STREAMS  ALL 

SELECT  ROADS  ALL 

ERASE  -  FLOOD  #  COLOR  8 

SHADE  # 

plot  #  #  #  #  color  14  2  1 

Stepg 

desc  coal. samp  at 

4 

1 

BSEARCH  COAL. SAMP 

SULFUR  GT  .5 

ERASE  -  PLOT  # 

LEGEND  #  LABEL 

3 

SULFUR 

Steph 

SAVE  # 

MYCOAL.PTS 

SELECT  MYCOAL.PTS  AL 

COMPUTE  # 

ELEVATION  -  TOT.OVERB 

COAL. TOP 

REPORT  # 

C0AL.RPT2 

YES 

NO 

3 

ELEVATION 

3 

TOT.OVERB    __ 

3 

1 

SYSTEM  SLIST  C0AL.RPT2 
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Answers:  Coal  Exercise  2 

Step  A 


SELECT  MYCOAL.PTS  ALL 
SELECT  S J. LEASE  ALL 

wiNDCfW  #  #  [4t  =  Active  table  ID  number] 

PLOT  #  # 

ZOCai  POINT 

.8 

/*  crosshairs 

ERASE 

PLOT  #  # 

GENERATE 

S J. BORDER 

MYCOAL.PTS 

EXERCISE2 

1990 

BODNDARY  FOR  COAL  HOLES 

YES 

1 

3 

C 

R 

1 

/*  crosshairs 

Q 

Stepb 

SELECT  S J. BORDER  ALL 
WINDOW  # 
ERASE 
PLOT  #   # 
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Answers:  Coal  Exercise  3 

Step  A 


FREE  ALL 

SELECT  S J. BORDER  ALL 

WINDOW  #  [#  =  Active  table  ID  number] 

SELECT  MYCOAL.PTS  AT 

3 

COAL. THICK 

1 

VARIOGRAM  # 

3 

COAL. THICK 

TITLE 

NO 

NO 

NO 

TES 

1000 

5 

22 

0 

22 

22 

22 

45 

22 

60 

22 

90 

NO 

YES 

coal.vario 
Stepb 

TESTGRID 

M 

250 

250 

N 

PLOT  # 
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Stepc 

bsearch  mycoal.pts 
coal. thick  gt  0 

GRID  # 

GRD.8PT 

3 

/*  2  carriage  returns  for  defaults 

1 

GRID  FOR  SJCOAL 

M 

250 

250 

N 

3 

COAL. THICK 

22 

GRID  # 

GRD.QS 

5 

/*  carriage  return  for  default 

1 

GRID  OF  COAL.PTS  QS 

M 

250 

250 

N 

3 

COAL. THICK 

7 

N 

DESCRIBE  GRD.8PT  HEADER 

DESCRIBE  GRD.QS  HEADER 

MAPS 

RENUMBER  GRD.QS  FOR  GRDR.QS  AS  3  TO  -9999  TH  3  AS  55  TO  55  TH 

9999 

BTB 

SELECT  GRD.8PT  ALL 

SELECT  GRD.QS  ALL 

PLOT  #  #  COLOR  2  4 

5,0,55 

5,0,55 
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Answers:  Coal  Exercise  4 

Step  A 


SELECT  GRD.8PT  ALL 

WINDOW  #[#  -  Active  table  ID  number] 

SELECT  MYCOAL.PTS  ALL 

ERASE 

PLOT  # 

LEGEND  #  LABEL 

3 

COAL. THICK 

ERASE 

PLOT  #  COLOR  5 

2,4,54 

PROFILE  #  /*  use  crosshairs  to  create  profile 

ERASE 

MAPS 

3D  GRD.8PT  RO  45  AB   30   INT   1  MAG  2 

ERASE 

BYE 
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Answers:  Coal  Exercise  5 

Step  A 


SELECT  MYCOAL.PTS   ITEM 
\1 

SAVE  #  [#  =  Active  table  ID  nnmbeT] 

MYPTS.NEW 

SELECT  MYPTS.NEW  AT 

3 

COAL. THICK 

1 

GRID  # 

SJTHK.GEU3 

3 

/*  two  carriage  retiims  for  defaults 

1 

GRID  FOR  SJCOAL 

M 

250 

250 

N 

3 

COAL. THICK 

22 

Steps 

describe  sjthk.grd  header 

MAPS 

RENUMBER  SJTHK.GRD  FOR  SJTHK.GRDR  AS  0  TO  -999  TH  0 

Stepc 

DELETE  SJTHK.GRD 

RENAME  SJTHK.GRDR  TO  SJTHK.GRD 

BYE 

Stepd 

SELECT  SJTHK.GRD  ALL 

SELECT  S J. LEASE  ALL 

ERASE 

PLOT  #  #  #  COLOR  5  4  2 

2,4,54 
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Answers:  Coal  Exercise  6 

Step  A 


SELECT  SJTHK.GRD  ALL 

WINDOW  #[#  =  Active  table  ID  number] 

SELECT  S J. LEASE  ALL 

POLYCELL  # 

PC. LEASE 

POLYCELLED  LEASE 

M 

250 

250 

2 

2 

DESC  SJTHK.GRD  HEADER  -  DESC  PC. LEASE  HEADER 

MAPS 

WINDOW  SJTHK.GRD 

PAGE  3000 

TOTAL  SJTHK.GRD  BY  PC. LEASE  FACTOR  1  FOR  LEASE. RPT 

BYE 

SYSTEM  SLIST  LEASE. RPT 

AREA  # 

/*  calculate  tonnage  using  vector  acreage 

/*  (area)  *  (avg.  thickness)  *  (factor)  =  tonnage 

/*  (9933.02)  *  (22.551)  *  (1770)  =  396,479,175  tons 

Stepb 

select  landlines  all 

OVERLAY  #  # 

OV.  SECTION 

YES 

0 

16 

1 

OVERLAY 

1991 

SJCOAL 

OVERLAY  OF  LEASE  AND  SECTIONS 

YES 

SELECT  OV. SECTION  ALL 

POLYCELL  # 

PC. SECTION 

SECTIONS  WITHIN  THE  LEASE      .  . 
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M 

250 

250 

2 

2 

MAPS 

TOTAL  SJTHK.GRD  BY  PC. SECTION  FACTOR  1  FOR  SECTION. RPT 

BYE 

SYSTEM  SLIST  SECTION. RPT 


OPTIONAL 


AREA  #  [AREA  the  vector  map  OV.  SECTION] 

/*  use  calculator  to  determine  tonnage  per  section  within  the  lease,  using 
vector  acreage 


SUBJECTS  VECTOR  AREA  *  FACTOR        *        AVG.THK        -  TONNAGE 


01;T23N.R13W 

632.01 

1770 

21.881 

02;T23N.R13W 

627.20 

1770 

24.854 

04;T23N.R12W 

685.38 

1770 

31.256 

05;T23N.R12W 

668.10 

1770 

43.716 

06;T23N.R12W 

665.03 

1770 

45.291 

07;T23N.R12W 

688.12 

1770 

30.247 

08;T23N.R12W 

702.85 

1770 

24.426 

09;T23N.R12W 

683.35 

1770 

21.514 

11;T23N.R13W 

662.08 

1770 

17.156 

12;  T23N.R13W 

658.02 

1770 

14.428 

13;T23N.R13W 

633.74 

1770 

15.218 

14;  T23N.R13W 

649.28 

1770 

14.227 

15;T23N.R13W 

643.10 

1770 

15.985 

17;T23VR12W 

671.31 

1770 

9.762 

18;T23"    R12W 

663.45 

1770 

14.529 

TOTALS  9933.01 
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Answers:  Coal  Exercise  7 

Step  A 


CCMPDTE   # 

TOT. OVERS  +  TOT.INTERB 

OVERINTB 

SELECT  MYCOAL.PTS  AT 

3 

OVERINTB 

1 

GRID  # 

OVERS. GRD 

3 

/*  three  ceurrlage  retxims 

GRID,  8PT,  FOR  OVERSDRDEN 

M 

250 

250 

N 

3 

OVERINTB 

22 


/#  =  Active  table  ID  nrnnbei] 


Stepb 

MAPS 

WINDOW  SJTHK.GRD 
MATH  {  OVERS. GRD  * 
(  SJTHK.GRD  *  .85 


(  1613.33  /  1770 
)  FOR  RATIO. GRD 


)    )    / 


Stepc 

2RASE 

FLOOD  RATIO. GRD  AS2TO0TH5AS3TO   5.00001  TH  10 

AS   4   TO   10.00001  TH   15  AS    11  TO   15.00001  TH   9999 

PLOT  PC. LEASE 

SHADE  PC. LEASE  AS  1  TO  1 

BYE 

SELECT  RATIO. GRD  AL 

PLOT  #   C   4 

5,0,15 
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Stepd 

MAPS 

RENDMBER  OVERS. GRD  FOR  FACTOR. GRD  AS  .85  TO  0  TH  100  , 
AS  .5  TO  100.00001  TH  9999 

MATH  {  OVERS. GRD  *  (  1613.33  /  1770  )  )  /  , 
(  SJTHK.GRD  *  FACTOR. GRD  }  FOR  RATI0GRD.2 
ERASE 

flood  rati0grd.2  as2to0th5as3to  5.00001  th  10  , 
as  4  to  10.00001  th  15  as  11  to  15.00001  th  9999 
plot  pc. lease 
shade  pc. lease  as  1  to  1 

Stepe 

RENUMBER  RATIO. GRD  FOR  COALPOT.GRD  AS  1  TO  0  TH  5  , 

AS  2  TO  5.00001  TH  15  AS  3  TO  15.00001  TH  9999 

ERASE 

FLOOD  COALPOT.GRD  AS2T01AS3T02AS4T03 

PLOT  PC. LEASE 

SHADE  PC. LEASE  AS  1  TO  1 
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Answers:  Coal  Exercise  8 

Step  A 


SELECT  SJTHK.GRD  AL 

WINDOW  #  [#  =  Active  table  ID  number] 

SELECT  COAL. SAMP  AT 

3 

SULFUR 

1 

SELECT  COAL. SAMP  AT 

3 

ASH 

1 

GRID  # 

SULFUR. GRD 

5 

<CR> 

1 

GRID  FOR  SJCOAL  SULFUR  PERCENTAGE 

M 

250 

250 

N 

3 

SULFUR 

5 

N 

GRID  # 

ASHGRD 

5 

<CR> 

1 

GRID  FOR  SJCOAL  ASH  CONTENT 

M 

250 

250 

N 

3 

ASH 

5 

N 

DESCRIBE  SULFUR. GRD  H  -  DESCRIBE  ASHGRD  H 

MAPS 

WINDOW  SJTHK.GRD 


AppMMilx       87 


RENUMBER  ASHGRD  FOR  ASH.GRD  AS  0  TO  -50  TH  0 

DELETE  ASHGRD 

CROSS  SULFDR.GRD  WITH  ASH.GRD  FOR  MINABIiE  , 

AS  1  TO  0  TH  .6  AND  0  TH  20  , 

AS  0  TO  .60001  TH  9999  AND  20.00001  TH  9999 

FLOOD  MINABLE  AS  4  TO  1 

ERASE 

Stepb 

cross  minable  with  ratio. grd  for  smcoal  , 
as  1  to  1  th  1  and  0.5  th  5 
flood  minable  as  4  to  1 
shade  smcoal  as  3  to  1 

ERASE 

Stepc 

CROSS  PC. LEASE  WITH  SMCOAL  FOR  SMCOAL. LSE  , 

AS  1  TO  1  TH  1  AND  1  TH  1 

LABEL  SMCOAL. LSE  SU 

16 

3 

1 

<CR> 

SMC  IN  LEASE 

<CR> 

1 

FLOOD  MINABLE  AS  4  TO  1 

SHADE  SMCOAL  AS  3  TO  1 

SHADE  SMCOAL. LSE  AS  9  TO  1 

COLOR  2 

PLOT  PC. LEASE 

COLOR  1 

CONTOUR  SJTHK.GERD  INT  1 

S7  £P  D 

TOTAL  SJTHK.GRD  BY  SMCOAL. LSE   LABEL  TONS   FACTOR   1 
TOTAL  SJTHK.GRD  BY  SMCOAL. LSE   LABEL  TONS   FACTOR   17'  3 
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Answers:  Coal  Exercise  9 

Step  A 


describe  dhdem.m  header 
describe  dhdem.m  projection 

Step  b  through  f 

MAPS 

WINDOW  DHTHK.GRD 

CUT  DHDEM.M  WINDOW  FOR  DEM. CUT 

RESAKPLE  DEM. CUT  FOR  DEM. RES  SIZE  250 

MATH  INT  (  DEM. RES  *  3.28083989  )  FOR  DEM. FT 

DESCRIBE  DEM. FT  HEADER 

DESCRIBE  DHTHK.GRD  HEADER 

LABEL  DEM. FT  HEADER 

DEM  IN  FEET 

OZ 

1991 

WILSON 

1991 

Stepg 

DELETE  DEM. CUT  ~ 

DELETE  DEM. RES 

BYE 

BYE 
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Answers:  Coal  Exercise  10 

Step  A 


SELECT  DEM. FT  ALL 

WINDCW  #[#  =  Active  table  ID  number] 

SELECT  DH.PTS  ALL 

POLYCELL  # 

PC. COLLAR 

POLYCELLED  COLLAR  ELEVATION 

M 

250 

250 

2 

6 

3 

COLLAR. ELV 

POLYCELL  # 

PC.SDBJ 

POLYCELLED  POINT  MAP,  SUBJECT 

M 

250 

250 

2 

2 

Steps  b  and  c 

MAPS 

WINDOW  DEM. FT 

MATH  PC. COLLAR  -  DEM. FT  FOR  DIFF 

CATEGORIZE  DIFF  FOR  DIFF.T? 

PAGE  3000         — 

TOTAL  DIFF.T7  BY  PC.SDBJ  LABEL  DIFF.RPT 

SYSTEM  SLIST  DIFF.RPT 

Stefd 

delete  pc. collar 

DELETE  DIFF.T7 

DELETE  PC.SDBJ 

DELETE  DIFF.T7M 

DELETE  DIFF 

DELETE  DIFF.T7P 

BYE 

FREE   ######TH# 
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ft  « 


OPTIONAL 

MAPS 

CATEGORIZE  DIFF  FROM  50  TH  9999  FOR  DIFF.T7P 

CATEGORIZE  DIFF  FROM  -9999  TH  -50  FOR  DIFF.T7M 

PAGE  3000 

TOTAL  DIFF.T7P  BY  PC.SUBJ  FOR  PLUS.RPT 

TOTAL  DIFP.T7M  BY  PC.SUBJ  FOR  MINUS. KPT 

SYSTEM  SLIST  PLUS.RPT 

SYSTEM  SLIST  MINUS. RPT 

BYE 

BYE 
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Answers:  Coal  Exercise  1 1 

Step  A 


MAPS 

WIKDOW  DHTHK.GRD 

MATH  DEM. FT  -  DHSTR.GRD  FOR  OVERB.DEM 

DESCRIBE  OVERB.DEM  HEADER 

RENUMBER  OVERB.DEM  FOR  0VERB.DEM2  AS  0  TO  -9999  TH  0 

FLOOD  OVERB.DEM  AS  7  TO  -9999  TH  0 

CONTOUR  OVERS. DEM2  INT  100 


.<s^s--^^' 
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